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Abstract 

The properties of excited nuclear matter and the quest for a phase 
transition which is expected to exist in this system are the subject 
of intensive investigations. High energy nuclear collisions between fi- 
nite nuclei which lead to matter fragmentation are used to investigate 
these properties. The present report covers effective work done on the 
subject over the two last decades. The analysis of experimental data 
is confronted with two major problems, the setting up of thermody- 
namic equilibrium in a time-dependent fragmentation process and the 
finite size of nuclei. The present status concerning the first point is 
presented. Simple classical models of disordered systems are derived 
starting with the generic bond percolation approach. These lattice 
and cellular equilibrium models, like percolation approaches, describe 
successfully experimental fragment multiplicity distributions. They 
also show the properties of systems which undergo a thermodynamic 
phase transition. Physical observables which are devised to show the 
existence and to fix the order of critical behaviour are presented. Ap- 
plications to the models are shown. Thermodynamic properties of 
finite systems undergoing critical behaviour are advantageously de- 
scribed in the framework of the microcanonical ensemble. Applica- 
tions to the designed models and to experimental data are presented 
and analysed. Perspectives of further developments of the field are 
suggested. 
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1 Introduction 

Protons and neutrons are able to form bound nuclei with well defined 
lifetimes. They are finite and contain a relatively small number of nucleons 
which varies from one to several hundreds. The nucleons experience a strong 
short range interaction which acts together with the long range Coulomb 
interaction between the charged protons. The interactions between nucleons 
determine the ground state properties and the low energy excitation spectrum 
of nuclei which have been studied over several decades. Their description by 
means of very different collective and microscopic, semi-classical and quan- 
tum models has lead to a satisfactory understanding of their properties [T|-|3| . 

With increasing excitation energy the detailed information coded in the 
spectrum and the wave functions of the system gets less and less meaningful 
and stringent for a satisfactory description of the system. Statistical argu- 
ments get relevant. Random matrix theories, for instance, have proved to be 
a very efficient tool in this regime |4],|[]. Under certain circumstances experi- 
ments reveal that the behaviour of nuclei can be described in the framework of 
thermodynamics. Excited nuclei behave like equilibrated compounds which 
may decay through emission of gamma rays, nucleons or light clusters at 
energies corresponding to unbound states of the system ||. The behaviour 
of nuclei in this regime and its description in terms of thermodynamic and 
statistical concepts leads to intriguing questions related to the general prop- 
erties of thermodynamically equilibrated nuclear matter, either finite like 
nuclei or (quasi) infinite as it could be found in neutron stars, in the absence 
of the gravitational interaction. These questions are also triggered by the fact 
that nucleons are fermions and nuclear matter in its ground state, at zero 
temperature, shows the properties of a Fermi liquid of strongly interacting 
particles |7|-[10| . For high enough excitation energy it seems sensible to think 



that such a quantum liquid may go over to a quantum gas by undergoing 
a phase transition which would be observable in a infinite system and for 
which one would find characteristic signals in finite nuclei. 

This question has become a major source of interest over the last two 
decades in the nuclear physics community. The only way to get an experi- 
mental answer to it goes through energetic collisions between nuclei or nu- 
cleons and nuclei with the aim to generate excited finite nuclear matter. The 
analysis of the experiments must then be able to reveal its properties and 
eventually deliver signs for the existence of one or several phase transitions 
in the infinite system. 



Such ambitious objectives have led to intense experimental and theoreti- 
cal investigations. Several reports have been written on the subject over the 
last ten years. In the present review we aim to discuss two aspects concern- 
ing excited nuclear matter. The first concerns the conditions under which 
it is experimentally generated. The second deals with specific microscopic 
models which we think allow a realistic study and description of experimen- 
tal facts under the prerequisite that the physical conditions imposed by the 
experiment match those which are imposed by the theory. 

The content of the present review is the following. In section 2 we describe 
the status of the subject at the beginning of the early 80's. We present and 
discuss the attempts to interpret experimental results related to energetic 
heavy ion collisions and the theoretical models which were aimed to describe 
the outcome of such collisions, the fragmentation of the involved nuclei. We 
then turn over in section 3 to the two major points which have to be taken in 
account in order to be able to compare experimental results with the type of 
models which will be introduced. These points concern the finiteness of the 
system and the problems related to its properties with respect to the concept 
of thermodynamic equilibrium. Section 4 is devoted to the outcome of very 
energetic nuclear collisions. Nuclei can break up into pieces. Fragments with 
different sizes, mass and charge numbers, are generated. With the presently 
available detection instruments they can be recorded. Their charge and, to 
some extent, their mass can be measured. Charge and mass distributions as 
well as related quantities are a priori non-trivial observables which contain 
interesting information about the fragmentation process and the behaviour 
of excited nuclear matter. We present and discuss the concepts and models 
which have been introduced in order to describe the fragmentation character- 
istics of nuclei, as well as the possible reasons for which rather simple-minded 
models describe these observables with a remarkable success. In section 5 
we introduce microscopic statistical models which may behave as generic ap- 
proaches for the study of the thermodynamic properties of excited systems 
composed of interacting nucleons. These models have either been borrowed 
from other fields of physics or consist of adapted derivatives of these models. 
They describe systems which show phase transitions and whose outcome can 
be confronted with experimental data. Since they are finite, tests which en- 
able to characterize the phase transitions have been proposed. Some of these 
tests are introduced in section 6 and their application both on theoretical 
models as well as experimental results are presented and discussed in the 
perspective of the existence of a phase transition in excited nuclear matter. 
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In section 7 we develop a summary, comment the present status of the field 
and suggest perspectives which may lead to further progress in the field. 

Nuclear fragmentation and the search of phase transitions is a very ac- 
tive field of research. This is best seen by looking at the number of preprints 
and publications which appear regularly. We tried to cover the newest re- 
sults concerning the aspects of the subject with which the present report is 
concerned, extending to summer 2000. We also tried to be as exhaustive as 
possible in the presentation and quotation of the work which has been done. 
We apologize for the contributions which could have been forgotten because 
of lack of information about their existence. 



2 Properties of excited nuclei and nuclear 
matter : first experimental and theoretical 
attempts 

Finite nuclei in their ground state and low-lying excited states have been 
studied for over 50 years. They are still objects of interest, in particular the 
so called exotic nuclei which contain a larger number of protons or neutrons 
than those which lie in the so called valley of stability These objects can 
be conceptually considered as finite representatives of infinite nuclear matter 
which may be realized in large objects of our universe like neutron stars where 
nucleons experience of course both the nuclear and gravitational interactions, 
the last one being absent from our present considerations. 

Extensive theoretical studies of the ground state properties of infinite 
nuclear matter have been pursued over several decades. They were essentially 
aimed to determine the non relativistic nucleon-nucleon interaction which 
originates from meson exchanges between nucleons. Its precise expression can 
be fixed by means of nucleon-nucleon scattering amplitudes jll-15] . Different 



expressions of this interaction have been introduced in order to construct 
effective interactions used in detailed spectroscopic studies of the ground 
state and low excited states of nuclei, see f.i. refs. |TB,[I7|. 



At high enough excitation energy, nuclei which form liquid drops of 
fermions get unstable. It is tempting to trace a parallel with the behaviour 
of macroscopic liquids and to conjecture that the nuclear liquid may go over 
into vapour and that there may exist a transition from a liquid to a gas un- 
der precise thermodynamic conditions, in the limit of an infinite system. But 
this is not necessarily the case, the evolution with increasing energy could 
well correspond to a smooth crossover from a bound to an unbound system 
of interacting particles. The question whether nuclear matter may exist in 
different phases is one of the important questions raised by the nuclear com- 
munity and a central point of the present review. The first serious attempts 
to answer this question started about two decades ago. 

2.1 Critical phenomena in nuclear matter : experi- 
mental facts and early interpretations 



2.1.1 Fragmentation of nuclei 

Experimental information about excited nuclear matter can only be gai- 
ned through accelerator induced nuclear reactions by means of which excited 
nuclei can be generated and studied. Energetic beams of particles or heavy 
ions shot on nuclei lead to their fragmentation into pieces. The final multi- 
plicity of fragments of different sizes which are collected by the detectors is 
correlated with the degree of violence of the collision. To our knowledge, the 
first experiment which was concerned with the possible observation of signs 
related to a phase transition in nuclei was performed in 1982 by Minich et 
al. | I5fl . Protons with energies between 80 and 350 GeV were shot on xenon 



and krypton targets. The resulting experimentally detected fragments with 
A nucleons (12 < A < 31) and Z protons led to fragment yields Y(A,Z) 
which could be parametrized in the form 

Y(A,Z)cxA-Tf(A,Z,T) (1) 

where r is a positive exponent and f(A, Z, T) a Boltzmann factor which de- 
pends on a temperature T and further constants which were fixed by a fit 
procedure to the experimental data. Expression ([]]) was in fact inspired by 
Fisher's droplet model |l9j which for the pupose of the analysis was general- 



ized to two types of particles, neutrons and protons. The value of the critical 
exponent extracted from the fit of the data was r = 2.64 and 2.65 for xenon 
and krypton respectively. It has to be compared to the value r = 2.33 ob- 
tained in mean field theory p0 |. Fisher's condensation theory predicts that 



the transition from a liquid to a gas proceeds via the formation of droplets 
whose size distribution follows a power law at the critical transition point. 
The authors do not claim that these results would be a proof for the existence 
of a critical behaviour but that they are consistent with it. Following the 
droplet model, the increase in excitation energy would generate instabilities 
in nuclear matter and would lead, at some critical temperature T c , to the 
disassembly of the system, which is here finite, into smaller pieces of all sizes 
which are experimentally detected. 

This pioneering experiment has been prolongated in further work which 
was pursued by Panagiotou et al. GTLE2]. The authors introduced a system- 



atic study of twelve different high energy reactions induced by protons and 
carbon ions on Ag, Kr, Xe and U targets. They measured fragment yields 
corresponding to nuclei with 3 < Z < 22. Using the same point of view as 
in ref. |[H|] the yields of fragments of size A were parametrized in the explicit 
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form 



where 



and 



P(A) oc A- k exp{[-a s (T)A 2 / 3 -a v (T)A + fi(T)A]/T} 

A-ky-A 2 l 3 Y A 



X = exp[-a s (T)/T] 



Y = exp[-(a„(T) - n{T))/T\ 

Here T is the temperature of the fragmenting system, a s ,a v ,fi are the surface, 
volume free energy per particle and the chemical potential respectively. T was 



determined independently for the different systems by different means [21 



Under the assumption that there exists a critical temperature T c at which 
a s (T) vanishes, the different quantities were further parametrized as 



a s (T) 
a v (T) -/x(T) 



18.4(1 - T/T c 
6(1-T/T C ) 2 



where b is a coefficient which has to be fixed. For T < T c where the gas and 
liquid coexist 

a v (T) - n(T) = 



at T = T r . Furthermore 



hence Y 



a s {T c 



hence X = 1 

For T > T c one assumes that a s (T) ~ 0. 

The expression P(A) was fitted to the experimental data in order to fix 
b and T c , for different values of k (=1.6 - 2.33) and in a temperature range 
T' — T ± 1 where T is the supposed common temperature of all emitted 
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fragments. The values of T were chosen within an interval of 2 MeV in 
order to take the uncertainty related to this quantity in account. If one fixes 
the critical temperature to the lowest value of t(T) in a parametrization 
P(A,T) oc A~ T ( T ^ where r is a so called apparent exponent, then the fitted 
value of T c coincides for 1.7 < k < 1.8 with T c ~ 12 MeV. Further analysis of 



the same type [22] confirmed these results, i.e. a critical temperature which 



lies in the range 9.5 - 13.5 MeV and an exponent k smaller than 2. 

2.1.2 Thermodynamic interpretation of fragment yields : com- 
ments and discussions 

The model which is aimed to reproduce the outcome of the fragmentation 
reactions presented above presupposes that the process can be described in 
the framework of thermodynamics and the grand canonical ensemble. It 
assumes that the fragmenting system is in thermodynamic equilibrium. This 
assumption is strong since it implies that there exists a relaxation time which 
is not longer than the time interval which separates the beginning of the 
process and the time at which fragments are formed. We shall come back to 
this point at different places below. 

The system undergoes a phase transition at a characteristic temperature 
T c . For T < T c a gas and a liquid coexist. The liquid phase is made of 
particles and composite droplets (fragments) whose energy is composed of 
a volume and a surface contribution. At T = T c the surface contribution 
disappears, the droplets disassemble into independent nucleons, hence form 
an homogeneous gas phase which is reached for T > T c . 

Taking the problem of thermal equilibrium apart, further thinking about 
this attractive picture leads to the following comments. First, the data anal- 
ysis is restricted to a finite range in the size of fragments. Light and heavy 
species are not taken in account. This may introduce a bias which is not 
under control. It has to be noticed that the exponents r and k which were 
obtained in refs. [2~I] and [2~2"j do not coincide with the exponent expected 



from the Fisher model which describes classical and infinite systems. Finite 
size effects should play a role in the determination of critical exponents, here 
they are not taken in account. Finally fragments and particles are free, i.e. 
they do not interact with each other. Even if the nuclear interaction is weak 
at the fragmentation stage because the system has already expanded, the 
Coulomb energy is present and certainly non negligible. Hence it should be 
taken into account. 
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At this stage, even though there were encouraging signs for the exis- 
tence of a critical behaviour of nuclear matter, it appeared that the observed 
phenomenon was not reducible to an ordinary liquid-gas phase transition as 
described by Fisher's droplet model. However, these first steps and the con- 
cepts they introduced paved the way to a huge amount of theoretical work 
which we shall briefly describe below, as far as the existence of different 
phases is concerned. 

2.2 Microscopic interpretations 

Here we restrict ourselves to those approaches which have been introduced 
as stationary and non-stationary descriptions of finite excited systems. They 
were aimed to describe the response and fate of nuclei when energy is im- 
parted to them. 

2.2.1 Temperature-dependent Hartree-Fock approaches 



The droplet model [19] is a classical and phenomenological description 
of an excited, inhomogeneous and unstable matter phase. The most com- 
mon approximation introduced in the framework of the quantum many-body 
problem is Hartree-Fock theory in which nucleons move in a self-consistent 
mean field which can be described by an effective density-dependent interac- 
tion. This framework has been used up to the present time for the study of 
static nuclear properties, a detailed review of its applications can be found 



f.i. in |23j and references therein. The theory has been extended to the 



study of systems at finite temperature in the framework of the canonical 
and grand canonical ensembles p4|-f29[1 . Using point contact and density- 



dependent effective interactions like Skyrme interactions it is easy to derive 
explicit expressions of different equations of state such as the expression of 
the pressure P as a function of the density p for different values of the temper- 
ature T. The typical behaviour of these quantities is presented in Fig. [I]. One 
observes the existence of an area limited by an envelope, the spinodal line, 
where the pressure decreases with increasing density. This is the signal for 
the presence of an instability of matter which corresponds to the formation of 
an inhomogeneous medium. It is interpreted as a zone where liquid and gas 
coexist and which is not thermodynamically accessible. The transition from 
an homogeneous to the unstable regime corresponds to a first order phase 
transition in the thermodynamic limit. The transition is of second order at 
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the critical point corresponding to the critical temperature T c , at the maxi- 
mum of the spinodal line where d 2 P/dp 2 = 0. This behaviour establishes a 
connection with the droplet model physics. Condensation has been studied 
in this framework |20|| . The calculations fix the critical density p c and the 



critical temperatures. The presence and importance of the Coulomb interac- 



tion was first taken into account by Levit and Bonche |50| who introduced 
a description in which the system is composed of a bound, excited nucleus 
surrounded by an external gas made of light particles in thermal equilibrium 
with the nucleus. The presence of the long range Coulomb force produces a 
sizable qualitative effect. It results in the apparition of a so called limiting 
temperature, T/ im , which restricts the instability zone to temperatures which 
are lower than T c and hence introduces an important modification. Similar 
investigations were performed in refs. |H|,|3l| with different Skyrme interac- 



tions and refinements concerning, in particular, the treatment of the vapour 
charge. 

The temperature-dependent Hartree-Fock approach with density-depen- 
dent effective interactions is a simple and elegant approximation to the many- 
body problem of finite excited systems. A priori it is possible to find the 
correspondence with central ingredients of the droplet model like the concept 
of surface in the liquid phase. It is however difficult to control the degree of 
realism of these microscopic quantum approaches. The fact that the process 
is described in the framework of the canonical or grand canonical ensemble 
when the fragmenting system is closed and has a fixed and small number of 
particles may lead to difficulties in the characterization of the transition. We 
shall come back to this point in the sequel. There are no means to obtain 
fragment yields which could be compared to the experiment since there are 
only two types of species, a bound nucleus and light particles which constitute 
respectively the liquid and the vapour phases. This is of course directly 
related to a major difficulty which comes from the mean field description. The 
phases are homogeneous. Many-body correlations induced by the neglected 
residual two-body interactions are absent, but their effect is very important 
at the considered excitation energies. This may also explain the high values 
of the temperature T c ~ 15 — 20 MeV where the system becomes unstable. 
Indeed, the interpretation of experiments indicates that these temperatures 
could be much lower, of the order of 4 — 6 MeV. 

Hence this kind of approach is at most indicative for the determination of 
the properties of finite excited nuclear matter in thermodynamic equilibrium. 
It worked however as an incentive to develop more realistic approaches which 
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take explicitly care of the many-body aspect of the problem. We describe 
some of them below. 

2.2.2 Time-dependent descriptions of nuclear fragmentation 

The general microscopic framework describing energetic collisions leading 
to the decay of the system into species of all sizes is the quantum many-body 
scattering theory. In this fundamental approach the scattering S'-matrix 
elements Sfi measure the overlap between any arbitrary many-body initial 
scattering state i with any arbitrary final state / |33], |4j]. A full fledged 
microscopic approach of this type is however hopelessly out of scope because 
of its mathematical intricancy. Many pragmatic ways have been devised in 
order to circumvent this problem. They all rely on non relativistic formalisms 
in which non nucleonic degrees of freedom are parametrized by means of 
nucleon-nucleon potentials. All approaches which have been used in a time- 
dependent framework contain a minimum of classical ingredients. Indeed, 
all of them fix initial conditions at an initial time and follow the evolution 
of the system up to a final time in contradistinction with scattering theory 
which specifies both the initial and the final state of the system. 

It is not necessarily clear that classical concepts and assumptions should 
work in the present physical context, except for the fact that high energy 
processes are correlated with short wavelengths. One may however notice 
that processes like deep-inelastic or fission reactions which occur at lower 
collision energies can be reasonably and surprisingly well reproduced in terms 
of classical (or semi-classical) equations of motion |]35], [36[ . 



An explicit time-dependent description generally called "dynamic" de- 
scription shows a priori a further advantage, since it allows to avoid the 
crucial point concerning thermodynamic equilibrium. Indeed, the dynamical 
description of the evolution of the system may lead or not to thermodynamic 
equilibrium which follows a transient period of time in a collision process. 

There exists essentially two classes of models of this type. The first 
one transcribes the many-body equations of motion into kinetic equations 
of motion for the one-particle phase-space distribution function f(f,p,t) at 
point f with momentum p. Each phase space cell {[r, f + df], [p,p + dp\} 
evolves in time through the transport equation 

?l + l.% + W.% = m (2) 

dt m df df dp 
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where V(r) is a one-body potential and /[/], the so called collision term con- 
tains in principle all the information relative to 2—, 3—, n— body correlations 
which are generated by the two-body (possibly also 3—, . . . ,n— body) po- 
tentials beyond the average field V. This quantity and /[/] are in principle 
correlated and fixed in the framework of the initial many-body equation of 
motion in such a way as to satisfy conservation laws. In practice however 
they are determined through phenomenological considerations which preserve 
more or less the quantal aspects of the initial problem. These approximations 
are particularly tricky for /[/] which has to be chosen in such a way that 
it can be expressed in terms of one-body distributions, in order to ensure 
the closed form of @ pTjjSSfl. The mean field potential V is generally fixed 
independently, from purely phenomenological formulations to self-consistent 
density-dependent expressions derived by means of static Hartree-Fock cal- 
culations. As a consequence, a large amount of models of the type of Eq. (Ef) 



have been proposed in the litterature, see |37| and refs. therein. They have 
been used as descriptions of many-body systems, among which high energy 
nuclear collisions. It is not our aim to give a detailed report on the achieve- 
ments of these descriptions. It has been done in different publications and 
reports, see the references above. Here we want to restrict the subject to 
the discussion of some aspects of nuclear fragmentation which have been ad- 
dressed in this framework. They are of interest for the understanding of the 
fragmentation process and problems related to thermodynamic equilibrium, 
the possible existence and (or) the coexistence of phases in nuclear matter. 
In the formulation of Eq. (|2|) the mean field V plays an essential role 
and its choice is central for a realistic description of the evolution of the 
system at high energy. As we have seen in section 2.2.1, realistic mean field 
calculations in the framework of Hartree-Fock theory with effective density- 
dependent potentials V and the equation of state P(p) lead to a zone of 
instability and possibly of metastability whose separation from pure "liquid" 
and "gas" phases corresponds to the existence of a first order phase transi- 
tion, see Fig. [1|. Transport equations may be considered as time-dependent 
extensions of mean field approaches which furthermore avoid the problem of 
the existence or not of thermal equilibrium in the expanding and decaying 
system. Much work has been done in this framework in order to identify the 
existence and consequences of an instability zone where the initially homo- 
geneous system gets inhomogeneous and appears as formed of clusters which 
can be identified [^, ^0[. This has indeed been achieved in a pure mean 



field approximation, in the absence of the collision term, when the mean 

16 



field is supposed to possess a fluctuating contribution |l I -fIJJ , and (or) when 
the collision term induces fluctuations in the one-body distribution function 



f(r,p, t) [fiE| , |l?|j . These fluctuations develop collective modes like zero sound, 



some of which are unstable and break the system |]48|,|49||. The expansion of 
the system affects also the development of the instabilities. The presence of 
two types of particles, protons and neutrons has been taken into account [5C| 



and comparisons with experimental results have been made ||51|| . Quantum 
effects on the expansion dynamics have been investigated recently. It comes 
out that they are non negligible, especially when the energy of the collective 
modes which develop in the system are larger than the temperature p2| , [53 



This approach which describes the fragmentation of nuclei as a mechanical 
process shows however limitations due to the already mentioned fact that 
Eq. (H) is an approximation to the many-body problem whose expression 
allows for many more or less controlled phenomenological ingredients, such as 
stochastic contributions in the mean field and (or) the collision term. These 
are supposed to simulate the missing many-body features. Unfortunately it 
is difficult to control their degree of consistency, much physical information 
is hidden in numerics and may obscure the conclusions which can be drawn 
from comparisons with the experiment. On the other hand one can believe 
that transport models of the type given by Eq. (|2|) are able to describe the 
time evolution of global observables like inclusive reaction cross sections and 



collective properties like flow [37 



The second approach concerns molecular dynamics. It is also an ap- 
proximation to the time-dependent many-body problem. Different formu- 
lations starting from the most classical one, Classical Molecular Dynamics 
(CMD) try to pick up quantal aspects at different levels, going from "Quan- 
tum" Molecular Dynamics (QMD) which tries to approximate the many-body 
wave function at the initial time by means of coherent (gaussian) states and 
takes care of the Pauli principle |54|, [55| to the more recent Fermion Molec- 
ular Dynamics (FMD) [|56|-|58| and Antisymmetrized Molecular Dynamics 
(AMD) |59|-|62| which introduce the antisymmetrization of wave functions in 
a semi-classical way. 

In the simplest form (CMD) the time evolution of the system is described 
in terms of classical equations of motion for protons and neutrons which 
interact through a two-body potential representing the nuclear interaction 
and the Coulomb force which acts between protons. These interactions are 
implemented in such a way that the fundamental ground state properties of 
real nuclei are reproduced. The effects of the Pauli principle can be mim- 
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icked by means of a momentum-dependent potential [53]. Applications of 



this formalism have been developed into different directions. We shall come 
back in the forthcoming section to the description of heavy ion collisions in 
which theoretically calculated and experimental determined observables are 
confronted in order to fix the possible relationship between fragmentation 
and phases in nuclear matter. Here we concentrate on applications in which 
Molecular Dynamics (MD) has been implemented in order to study the prop- 
erties of a thermalized system for which the velocity distributions of particles 
are maxwellians. From there on the system evolves in time by expansion for 
different initial values of the temperature and the density f63|-|65| One follows 
the evolution of the temperature, the pressure and the fragment size distri- 
bution of the system with decreasing density. Fragments can be identified 
by means of different criteria which may be more or less realistic |6^, |66| , |67|| . 
In some cases the system crosses the spinodal line and enters the instability 
zone. Different observables show a behaviour which can be interpreted in 
terms of Fisher's droplet model, i.e. a fragment size distribution correspond- 
ing to a power law corrected for finite size effects and a corresponding critical 
exponent r ~ 2.2 to 2.3, which goes over to other characteristic shapes for 
other choices of the temperature as expected in the droplet theory, see Fig. |2|. 
These results along with the interpretation of other observables which will be 
defined below are interpreted as the outcome of a realistic simulation of the 
experimental situation and the sign for the existence of a critical behaviour 
of the system. The role and importance of the Coulomb interaction has been 
investigated [|64 |. Molecular Dynamics has also been used in order to con- 



struct the equation of state f65|, i.e. the thermodynamic properties which 



link density, pressure and temperature. It is again shown that there exists a 
critical point for finite systems, and a power law distribution with r ~ 2.23 
which can be interpreted in the framework of Fisher's model. 

However these results and conclusions concerning criticality have been 
put in question by another study |68| which claims that the observed signs 



have not to be interpreted as being related to a critical behaviour but for 
the existence of an intermediate regime. This behaviour would correspond 
to a trajectory in the (T, p) plane on which the evolving system does neither 
fall back to a drop nor develop into a gas. A power law distribution can be 
extracted from the data. The power law exponent depends sensitively on the 
initial temperature and the size of the system. In ref. |68|] it is claimed that 
the power law behaviour of mass yields expected for a system close to its 
critical point would be purely accidental. 
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Time-dependent transport descriptions of nuclear fragmentation seem a 
priori to be the most realistic way to the description of high energy nuclear 
collisions. They have been widely used in order to describe the possible out- 
of-equilibrium character of the fragmentation process which characterizes at 
least the first stage of the process |HJ. In the present section we concentrated 
on aspects related to a time-dependent description of thermodynamically 
equilibrated systems which expand and cool down, exploring different regions 
of the (P, p) and (T, p) planes. Calculations show that the results can be 
interpreted as the entrance of the system into a thermodynamic unstable 
regime. It seems however not clear how this can effectively be related to 
the existence of a phase transition which would characterize an equilibrated 
system heated up in a fixed volume. In an interpretation in which the system 
gets critical it is claimed that it is possible to extract critical exponents, at 



least the power law exponent r [531]. However this quantity is obtained from 
finite size systems and hence should be corrected for finite size effects. 

The time coordinate which is introduced in these descriptions raises the 
question of time scales, i.e. the question concerning equilibration in energetic 
light particle reactions or heavy ion collisions. This is the subject of section 
3. Last but not necessarily least, all these time-dependent approaches need 
a detailed description of the dynamics in terms of more or less sophisticated 
mean field and two-body potentials. It is always difficult to get a clear feeling 
for the sensivity of the outcome of numerical calculations and simulations 
to these details. Paradoxically, it is at criticality that the information on 
microscopic systems is the less sensitive to the details of the Hamiltonians 
which govern them. 

Finally it is worthwhile to mention here that it has been shown re- 
cently [[F(J that even far-from-equilibrium, spatially chaotic systems can show 
equilibrium properties such as ergodicity (equivalence between time and en- 
semble averages), detailed balance in microscopic processes, partition func- 
tions and renormalization group flow at coarse-grained lengths between mi- 
croscopic and macroscopic scales. This suggests that it might be possible 
to describe the global behaviour of some far-from-equilibrium systems in the 
framework of equilibrium statistical mechanisms. Whether this result could 
be some clue to the actual question raised in the present framework is an 
open and highly interesting question. 
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3 Thermodynamic equilibration and phase 
space descriptions of nuclear fragmentation 

If an isolated system of particles is such that there is no unbalanced force 
acting between its constituents and if it does not experience internal structure 
changes of its constituents it is said to be in mechanical and chemical equi- 
librium. Thermal equilibrium is realized if the relevant degrees of freedom 
which characterize the system share equally its macroscopic energy. One can 
then define a common temperature which characterizes the energy attributed 
to each degree of freedom. When all three conditions are realized the sys- 
tem is said to be in thermodynamic equilibrium Jn]] and can be described in 



terms of macroscopic coordinates, the thermodynamic coordinates like den- 
sity, pressure and temperature which do not change with time and are linked 
together through equations of state. In the present section we aim to present 
and discuss the difficult and for a large part unsettled problem of thermody- 
namic equilibrium in small systems like nuclei excited by means of energetic 
collisions between particles or nuclei called heavy ions if their mass is larger 
than 4. In section 3.1 we discuss questions related to the interval of time 
over which equilibrium can be reached in terms of model estimates. A con- 
frontation between experimental facts and phase space equilibrium models is 
presented in section 3.2. 

3.1 Thermodynamic equilibration in excited nuclear 
systems generated by nuclear collisions 

What is the scenario which governs energetic collisions? Data analysis 
and their interpretation show that several scenarii are at hand. Typically 
one expects that part of the bombarding energy is transferred to the internal 
degrees of freedom of the system, the system may fall into pieces, particles 
and fragments, and expand in time. The remnants are then collected at 
asymptotic distances in detectors which identify their charges and eventually 
their masses, along with their kinetic energies and distributions in space. 
The setting up of thermodynamic equilibrium over the whole or part of the 
interacting system of particles is a central but open question. Under favorable 
circumstances one may figure out that after some transient interval of time 
which starts at the beginning of the collision the system reaches at least 
partial equilibrium at some stage and then expands to infinity. 
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3.1.1 Low energy time scales 

At very low projectile energies the excitation of nuclei by means of light 
particles can lead to the formation of an excited complex whose behaviour 
can be interpreted within the framework of thermodynamics as a system in 
thermodynamic equilibrium and which can loose part of its excitation energy 
through the emission of particles, photons or through symmetric or asymmet- 
ric fission. The complex is the so called compound nucleus whose formation 
and decay properties were already studied in the 30's by N. Bohr [|T|.|7!2[ 
who conjectured that, in the framework of this scenario, the process can be 
explained as a two step mechanism in which the decay of the equilibrated 
system would be totally uncorrelated from the way it was generated. From 
a theoretical point of view compound nucleus theory is a phenomenologi- 
cal approach. It has been largely investigated and worked out in terms of 
microscopic models, among others random matrix theory. There exists an 
extended litterature on this subject which covers several decades of intensive 
work [73] . Indeed, the consequences of the model have been abundantly ver- 



ified by means of physical observables like excitation functions and angular 
distributions of emitted particles [|73|,[7j|. The lifetime r of the compound 



system can in principle be read from the compound nucleus decay width 
AE through AE ■ r > h leading to r ~ 10~ 20 — 10~ 17 s depending in par- 
ticular on the size of the system. These times are substantially larger than 
r ~ £/3 x lCT 23 s corresponding to the time necessary to light to cross a 
nucleus of diameter £ fm. 

At higher energies, above the Coulomb barrier which corresponds to bom- 
barding energies per particle E/A ~ 1 — 2 MeV, reactions involving heavy 
ions lead to a new mechanism the so called deep inelastic collision (DIC) 
process. Typical contact times during which the projectile and the target 
interact through the short range nuclear potential are much shorter than 
compound nucleus decay times, of the order of 10~ 21 s, not much larger than 
the contact time in so called direct reactions in which the ions strike each 
other peripherally in typical time intervals of 10~ 22 s. The interpretation of 
physical observables measured in these DICs shows that although the char- 
acteristic reaction times are quite short, the complex which is formed over 
the contact time of the two ions is sufficient to transfer sizable amounts of 
energy from the relative motion between ions into the excitation of the nucle- 
onic degrees of freedom and the characteristic features of the system can be 
interpreted in terms of an equilibrated system whose equilibration is reached 
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in as short as some units of 10 22 s f35|,[75|. 

3.1.2 Approaches to the equilibration problem at high energies 

How does this behaviour extrapolate to high energy events which can lead 
to the fragmentation of the system into many pieces of different sizes? The 
problem is intricate since it is difficult to guess the scenario (or the differ- 
ent scenarii) as already mentioned above. Hence one may proceed in steps, 
starting with simple-minded models. To our knowledge, the first attempt 
was made by Boal f76f , at the time where the experiments, which led to the 
belief that signs for the existence of a phase transition had been seen, were 
performed |H],^2|. In fact the question which was raised did not concern the 
equilibration but the fragment formation time. Former experiments J77J had 
shown that the ratio between inelastic proton reactions <j{p,p') and proton- 
neutron exchange reactions a(p, n), cr(p,p')/a(p,n), is of the order of 2 at 
100 MeV incident energy after correction for the Z/N ratio (Z, N = proton 
and neutron number). This is evidence for the fact that chemical equilibrium 
is not reached, hence strictly speaking thermodynamic equilibrium cannot be 
either, see above. The author of ref. J7|] postulates that a hot system gener- 
ated at high temperature cools down to a fraction of this temperature and, 
from thereon, experiences particle coalescence leading to the formation of 
clusters which come out as fragments. Defining Ni(t) as the number of clus- 
ters of size i which exist at time t normalized to the total number of particles 
At, the initial conditions are supposed to be given by 

iV\(£ = 0) = P 

Ni(t = 0) = 2 < i < A T 

where p is the matter density. The system starts at temperature T from 
an homogeneous distribution of particles and its evolution is given by the 
postulated rate equations 

^£^^W-£« (ife (3) 

i,j ■* i 

where the first term on the r.h.s. of (|3]) corresponds to the coalescence of 
species of size i and j to species of size k and the second term to the decay 
of species of size k to species of size i. The coalescence and break-up rates 
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are given by 

a = 4vr ( -^— j / dvv 3 a(v) exp(-/j,v 2 /2T) 

where \i is the reduced mass of the colliding species, v their relative velocity 
and a(v) the corresponding cross section. The solutions of the system of 
coupled equations ([$]) reaches the asymptotic regime after 4 x 10~ 23 s. The 
results are shown in Fig. |3] by the histogram for the case of a reaction p + Kr 
at energies between 80 and 350 GeV. Their global trend agrees quite nicely 
with the experimental points. The asymptotic time is in agreement with the 
c(p,p')/a(p, n) ratio and the estimated rate of cooling of the excited source, 
~ 10 _23 s. One should mention here that the energy is very high and that 
thermal equilibrium is postulated since one defines a temperature T which 
is furthermore kept constant over the duration time of the process. The 
expansion of the system in space is not taken into account. In a further work 
Boal and Goodman abandoned the coalescence scenario of an initially 



homogeneous and expanded system. They proposed a break-up mechanism 
related to the entrance of the system into an instability region of the equation 
of state as already discussed in section 2. The collision of energetic protons 
(E p = 300 MeV) with nuclei was simulated by means of a simplified transport 
equation of the Boltzmann-Uehling-Uhlenbeck (BUU) type. The authors 
emphasize the fact that thermal equilibrium, requested by a classical collision 
time between nucleons which is short compared to the expansion time of the 
nucleon gas, may not be established at a time where the system is already 
dilute, a regime which is reached after a very short characteristic time which is 
smaller than 10~ 22 s. This goes in the same direction as the non equilibrium 
a(p,p')/a(p,n) which was experimentally found. Notice however that the 
energies involved here are very large, non nucleonic degrees of freedom are 
not taken into account and the transport description is somewhat schematic. 
The different models are used in order to explain the increase of entropy of 
the system when it crosses the spinodal line and enters the instability zone. 
This entropy agrees with the quantity extracted from the experiment. 

Some aspects of fragmentation related to time scales have been discussed 
by Gross and coll. [75 ]. The reaction process goes through a transient out- 



of-equilibrium phase which may ultimately lead to the formation of an equi- 
librated excited system. In order to estimate the transient time the authors 
use a BUU-type equation f80H33| which is a mean field description without 
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collision term as described in section 2. They argue that this description 
may be valid at the beginning of the collision process, up to and not further 
than the time when the system may break into pieces. At energies of 50 to 
60 MeV-A this happens after 50—60 fm/c for a reaction 96 Mo on 96 Mo. Then 
the system expands and apparently reaches equilibrium at about 200 fm/c 
due to the dissipation of energy in the nucleonic degrees of freedom. How far 
the equilibrium is reached is however not firmly established. Once the dis- 
tance between the different species gets larger than 3 fm, they do no longer 
interact by means of the nuclear potential. The system has then reached 
the so called freeze-out stage in which it is in chemical equilibrium. A phase 
space description of fragments and particles which move under their mutual 
Coulomb interaction can then be introduced. The early time behaviour is 
confirmed by experimental results fSl ], 



Recently relativistic transport equations 55 have been used in order to 



simulate the heavy ion reaction Au on Au which has been studied experimen- 



tally |86],£7]]. The calculations show that after the collision the excited system 
separates into a participant region of highly excited matter and a spectator 
which reaches thermal equilibrium, approximately after 2 • 10 _22 s. Indeed, 
a local temperature can be defined through the knowledge of fragment ki- 
netic energies. It characterizes a system which is mechanically unstable and 
decays into fragments. The values of the temperatures obtained for differ- 
ent beam energies compare well with those which are extracted from the 



experiments [BRBU 



Other decay scenarii borrowed from low energy reaction mechanisms have 
been introduced [|9"0"H95l- Fragments are formed sequentially in time, by 



means of the successive binary decay of the excited species which already 
exist at a given time step. The process starts from an equilibrated excited 
system. In ref. |JIJ this mechanism is supplemented by a boost in time 
due to the early compression of the system at the beginning of the process. 
Most of these approaches rely on the Weisskopf formulation of the break-up 
rates [|96| of an excited bound cluster of particles. In ref. [94| the binary de- 



cay is interpreted in terms of fission barriers, using a formulation introduced 
by Swiatecki |9^]. The use of the Weisskopf expressions raises problems 



which shed doubt about their physical relevance at high energies. First, they 
get unrealistically small at energies corresponding to temperatures above 
4 — 5MeV f98j . Second [[FJJ, each decay event in this formulation should be 
independent of the preceding one. This is only possible if the characteristic 
emission time exceeds the characteristic time over which the Coulomb emis- 
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sion barrier changes for the following emission. This time is of the order of 
10 _21 s, much longer than the Weisskopf emission time and not compatible 
with characteristic fragment formation times as estimated from the exper- 



iments [99]. Hence, although confrontation of mass yields obtained in this 



framework with experiments may seem satisfactory it is difficult to believe 



this type of description [100] for highly excited systems. It is of course sen- 
sible to believe that fragmentation of nuclei proceeds through a sequence of 
break-ups, but it is difficult to guess a priori or to derive from some micro- 
scopic model the realistic expressions of the rates which govern them when 
the excitation energy is as high as several MeV/A. 

The evolution of the fragment emission time has been studied recently 
by means of two-fragment correlation measurements over a large range of 
excitation energies in reactions ir~ and p on 197 Au [ 101 . The fragments are 
emitted from a unique source which is considered to be in thermal equilib- 
rium. Confrontation of the data with results obtained by means of classical 
trajectory calculations allows to follow the emission time which decreases 
with increasing excitation energy. The result is interpreted as evidence for a 
cross-over from surface emission over long times as described by sequential 
decay and fast bulk emission which is characteristic of a break-up process 
due to mechanical instability. 

It has also been looked for further direct insight into the time scale 



problem. In ref. ||102| central-impact-parameter events from the reaction 
40 At + 51 V with incident energies between 35 and 85 MeV ■ A were anal- 
ysed for the shape of the emitted particle and fragment events in momentum 
space. It is expected that fragments which are emitted quasi-simultaneously 
show a close to isotropic emission in the centre of mass system, hence are 
spherical in shape, whereas emission times which are long as it may be in 
the case of sequential emission would correspond to the existence of an elon- 
gated emission axis, hence an ellipsoidal shape. Experimental events were 
confronted with numerical simulations. It comes out that for energies below 
35 MeV -A the fragmentation process is slow. For higher energies, the experi- 
mental momentum distributions lie in between the predictions corresponding 
to a slow and a quasi simultaneous break-up process. Even though the result 
does not give any quantitative clue to the problem, it shows at least that 
the characteristic fragmentation time decreases with increasing energy. A 



fragmentation experiment induced by m Ca + 40 Ca at 35 MeV- A [103 



con- 



fronted with different models related to quick and slow processes confirms 
the fact that ordinary binary sequential decay may not be the appropriate 
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fragmentation mechanism. 

Other characteristic time scales can be introduced such as the time for 
the occurence of the first break-up of the system tq, and the average length of 



the interval of time between successive break-ups tff ||104|| . These times can 
be estimated by means of relative angle correlation measurements between 
breaking species. If the break-up time is slow, Coulomb repulsion will hinder 
the emission at small relative angle and low relative velocity. Simulations of 
the correlation functions indicate that tff corresponds to prompt decay for 
excitation energies exceeding 4 MeV/A. Estimates of tff are shown in Fig. |j. 
A long time interval To has been observed for excitation energies lower than 
3 MeV/A in peripheral collision events. Even for high excitation energies of 
the order of 5MeV/A r ~ 100 — 150 fm/c. This leads to characteristic 
times of 300 — 400 fm/c for complete break-up which are rather large for 
heavy nuclei, close to the Fermi energy. Recently further emission times have 
been deduced from IMF-IMF correlation functions obtained from hadron 



induced multifragmentation events 101 . One observes an evolution from 



typical emission times r ~ 500 fm/c at excitation energies E* / A = 2MeV 
to r ~ 20 — 50 fm/c for E*/A = 5MeV and above. The observation of 
the transition is interpreted as a transition from surface to bulk-dominated 
emission. The characteristic times are somewhat shorter than those obtained 
in ref. ||104|,|105|] for excitation energies lower than 6MeV/A. 



Albeit the numerous and different types of investigations which have been 
made up to now the arguments concerning a quick thermodynamic equili- 
bration of the system and relying explicitly on time remain the subject of 
controversies since time is not directly under experimental control. These 
controversies are in particular raised by time-dependent descriptions of the 
collision process, in the framework of the Quantum Molecular Dynamics 
(QMD) models |f>3],|106],|i07|. These are classical molecular dynamics models 



which, to some extent, take care of the quantum nature of nucleons. In re- 
cent applications of this type of models to specific reactions which have been 
investigated experimentally in some details like 129 Xe on Sn at 50MeV-A 
incident energy, the confrontation between calculations and experimental re- 
sults shows good quantitative agreement, except for the mass yields corre- 
sponding to fragments with more than 30 nucleons and too small kinetic 



energies for fragments and particles ||108||. It is argued that the system of 



interacting particles is far from equilibrium, there is no sign for a thermally 
equilibrated behaviour of the species which are generated during the early 
stages of the process. In a further study |109|| QMD calculations are con- 
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fronted with 50 MeV-A Xe + Sn data which were also analysed in the frame- 
work of the Statistical Multifragmentation Model (SMM, see below) |110| . 
The approach assumes thermodynamic equilibrium at the freeze-out. Exper- 
imental data related to the multiplicities of light charged particles (LCP), 
the average kinetic energy of LCPs and intermediate mass fragments (IMF) 
are very close to those obtained by means of QMD and SMM simulations. 
This result concerning the average transverse kinetic energy of IMFs and the 
total kinetic energy of LCPs contradicts two arguments which were supposed 
to help to distinguish between a fragmentation in an equilibrated and non 



equilibrated system [[111 



The preceding discussion shows once more that the experimental data 
related to fragmentation processes are ambiguous as far as the production 
mechanism is concerned. It raises the question why the outcome of differ- 
ent descriptions should be so close to each other when the process is over, 
i.e. at the freeze-out. At the present stage where different models lead to 
similar results which agree more or less satisfactorily with the experiment it 
is tempting to conclude that all of them contain part of the truth, at least 
to a sufficient degree as to agree with the measured quantities. Since there 
exists no trustful method to settle the point, one has to collect a maximum 
amount of informations, as exclusive as possible, and confront them directly 
with theoretical descriptions. In the next subsection we intend to show that 
in fact most of the data can be understood in terms of systems which reach 
thermodynamic equilibrium at the freeze-out. 

As a final remark, one may wonder why the occupation of phase space 
leading to equilibrium goes so fast, seemingly even faster at high than at 
low excitation energy. In a quantum mechanical microscopic description, the 
decay width of a state |/> into a state \i> at the energy E is given by the 
golden rule expression 

T(E) oc 2rt\V ft \ 2 p f (E) 

where Vfi is the interaction matrix element between the initial and final state 
and pf the density of states at the energy of the final state. For increasing 
energy E, Pf(E) will increase exponentially and if Vfi does not decrease as 
fast as or faster than pf which is certainly the case, T(E) increases with E. 
If we associate a transition time r(E) oc %/Y{E) to this process, one sees 
that the transition time between different states of the system can decrease 
very fast and hence drastically accelerate the path to equilibration. 
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3.2 Thermodynamic equilibrium : confrontation of ex- 
perimental data with phase space models 

Excited nuclei in thermodynamic equilibrium are described by means 
of so called statistical models. They generally describe a situation where 
the system is at the freeze-out, the stage at which particles and fragments 
present in the system are located at relative distances d from each other 
which are larger than 2 — 3 fm, so that the nuclear interaction between 
them is negligible. Before we show applications of this type of models to the 
analysis of experimental results we present and discuss some of their aspects. 

3.2.1 The first phase space models 

To our knowledge, the first attempts to describe strongly excited systems 
of particles and bound clusters were made in the beginning of the 80 's. The 
concept of statistical multifragmentation was first introduced in ref . [pj| . In 



1981, Randrup and Koonin proposed a thermodynamic phase space model 
113|| . They defined the partition function Z of a classical system in the grand 



canonical ensemble in which the number of particles A, N neutrons and Z 
protons, as well as the energy, is fixed in the average. If Aq is the average 
number of particles, one defines 

uj = In Z/A = \^ uj at 

AT 

where T = 1/2(N — Z). The multiplicity of fragments characterized by A 
and T is given by 

/ An \ f27nnA\ 3 ^ 2 
Uat = \ T r ° X J ' 1 Rfr 2 ) ' Zat ' ex P[~^(^ - \iA- vT)\ 

In the first term which is an effective volume, \ is a parameter of order 
unity which can be fixed by comparison with experiment, r (= 1.15 fm) 
the nuclear radius constant (R sys t — t A 1//3 ). In the second term which is 
related to the kinetic energy contribution m is the nucleon mass and /3 the 
inverse temperature. The third factor takes care of the internal energy of the 
fragments 



^2g%exp(-/3e 



(i) N 
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where g^ T = 2j^ T + 1 is the degeneracy of the excited level i and e^ T the 
corresponding energy. In the fourth term Vat is the gound state mass ex- 



cess which is taken from a liquid drop formula ||114|| , pi and v are Lagrange 
multipliers which are fixed in such a way that the total energy E , number 
of particles A and isotopic composition are fixed, 

A ° = w^ nZ 

where Z = / esxp[—j3(Ef — \iAf — vTf)] 

f 

and the sum over / extends over all final states which are characterized by the 
number of fragments, their mass numbers, isospin projections, internal exci- 
tations, positions and momenta. In this description, particles and fragments 
do not interact, hence their spatial location is not relevant. 

The first confrontation with experiment was made with an extension of 
this model in which unstable species decaying by particle emission were in- 
cluded |115|| . The fragmentation was supposed to consist of a quick explosion 



followed by a slow evaporation from unstable excited fragments. The light 
particle generation ratios d/p, a/p, 3 He/p, t/p as well as the pion rates tt~ /p 
and 7r + jp were worked out and showed good agreement with experiment. 

Nuclei are finite objects and closed systems, with a fixed number of par- 
ticles and a fixed energy. An approximate microcanonical description was 
introduced later | 116 | along with an extension which allows for the existence 



of different fragmentation sources generated in the collision process, which 
either "explode" for a large enough energy content or decay sequentially for 
low excitation energy. In this form the model has been used for explicit 



numerical simulations 117 



In a further step Koonin and Randrup [ 118|| implemented the model with 



a more rigorous canonical and microcanonical description in which the inter- 
action between fragments and particles is taken into account. The essential 
quantity in the microcanonical framework is the density of states for a system 
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of Ap particles with energy Ep and located in a volume Q 

p{Q, A, E) = J2 6 ( a f ~ A)5{E F - E) 

F 

where F is the set of fragmentation configurations defined by a set of variables 
{A n ,f n ,p n ,e n } (n = 1, • • • ,Np) which specifies the mass, charge, position, 
momentum and internal energy of the species labelled by n. The total energy 
is written as 



E, 



n=l 



Pr. 



2mA r 



B n + e, 



+ 7) 2_^ V nn> 



n'^n 



where B n is the ground state binding energy and 

e Z n Z n i 



Vnn' 



\ I n. ' n. f \ 



+ V nuc (\r n -r n t\) 



the potentials which act between the fragments and the particles. Appli- 
cations of the model which includes metastable states showed the effect of 
the finiteness of the system which is contained in a fixed volume and the 
importance of the interaction energy between the constituents. 

3.2.2 The Berlin and Copenhagen phase space models 

The first statistical theory on fragmentation of finite nuclei which in- 
cluded an exact treatment of the Coulomb energy between charged particles 
and fragments was developed in ref. |119||. The most popular and success- 



ful models are the Berlin model called MMMC (Microcanonical Metropolis 
Monte Carlo) p|,|T20HT22| and the already quoted SMM (Statistical Multi- 
fragmentation Model) CT , |T23| - [T25[| . A description closely related to MMMC 



in its spirit was developed in Beijing ||126 



Except for technical though important points concerning Monte Carlo 
event samplings which are used in the practical implementation of these 
models, the main difference between MMMC and the Koonin-Randrup model 
concerns the fact that the system which is at the freeze-out density does not 
experience evaporation of light particles in its final stage. The volume of 
the system is kept fixed and spherical, the Coulomb interaction between the 
constituents is treated rigorously. The model has been extensively used in 
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order to confront calculated and experimental observables such as fragment 
correlation functions at high energy ||127|| and also low energy data ||128|| . It 



has been tested for the possible existence of a first order phase transition, see 
below, section 5. Recently it has been extended to non spherical shapes of the 
volume and the consequences of this generalization have been studied [ |129 



The Copenhagen approach has been initially worked out in the frame- 
work of the canonical ensemble [|123| , |124| and later on a microcanonical for- 
mulation was developed ||110fl . The approach is rather close to MMMC, see 



ref . ||130|| . The main differences concern the fact that particles can evaporate 
from excited fragments, the Coulomb interaction is treated in the Wigner- 
Seitz approximation and the volume in which the system is enclosed at freeze- 
out is not kept fixed. The problem of equilibration cannot be raised in the 
framework of the model. The question has been investigated by means of 



simple-minded transport equations ||110|. The authors rely on the obser- 



vation of collective flow of matter which is generated and whose transverse 
energy carried perpendicularly to the incident beam axis by the particles and 
the fragments can be measured. In this framework the confrontation with 
experiment shows good agreement. The interpretation of the data indicates 
that fragments are generated very early in the process, over time intervals as 
small as 50fm/c. This shows that fragment formation goes as a break-up 
and not a condensation process if the description is realistic. The instability 
growth is accelerated by the collective expansion of matter when the system 
enters the spinodal region. It is argued that the energy transferred to the 
system is converted into heat, generating a chaotic system which thermalizes 
somewhat before break-up, hence after a short interval of time. 

3.2.3 Thermal equilibrium : confrontation of models with exper- 
imental facts 

Since there exists no direct way to measure the evolution and to find out 
experimentally whether thermodynamic equilibrium is reached or not, one 
has to rely on the confrontation between measured and model-dependent 
physical observables extracted from numerical simulations. This procedure 
has been followed a number of times in the recent past ||131| — |137f| . We give 



here some examples which show how far equilibration is at least compatible 

with the experimental outcome even if it is not possible to prove it rigorously. 

In an experiment m Ar + 58 Ni at different energies (52, 74, 84 and 

95 MeV ■ A) Borderie et al. |138| . |13Sj| analysed the outcome of fragmenta- 
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tion events in which 90 % of the nucleons are identified with the detector 
INDRA. They showed that these events correspond essentially to the va- 
porization of the system into light particles and clusters. The forward and 
backward spectra of particles are superimposable, which at low energy would 
be interpreted as a sign for the formation of a compound system. The en- 
ergy distribution functions of the different light species show experimental 
tails whose slopes are comparable to within 30 %. The authors introduce 
simple thermodynamic models like the so called Quantum Statistical Model 
(QSM) ||140[ - |142| which works in the framework of the grand canonical en- 



semble, in an attempt to reproduce the measured observables. A detailed 
comparison shows that all measured quantities are well reproduced, in par- 
ticular light particle yields, their average kinetic energies and the variances 
of multiplicity distributions. 

The fragmentation events observed in the reaction 129 Xe+ 197 Au at 30, 40, 



50 and 60MeV-A projectile energy have been analysed in some detail |143 
The average transverse energy (E t = £\ E,- t sin 2 6i , where E^ sin 2 Q^ is the 
projection of the energy associated with species i in the direction perpendic- 
ular to the incident beam) measures the energy transferred from the relative 
motion of the ions to the excitation energy of the system. The analysis shows 
that the average number of light charged particles < Nlcp > and the av- 
erage number of IMFs and their corresponding transverse energies increase 
monotonously with increasing E t as shown in Fig. [5| This contradicts the 



results of | |144| | where the transverse LCP energy saturates. This saturation 



can in fact be understood as being due to instrumental problems. The fea- 
tures observed in Fig. || can be interpreted and understood in the framework 
of a statistical decay mechanism ||134j| . The results can also be reproduced 



in the framework of phase space models like SMM ||110|| . 



Further analysis of the reaction 129 Xe + Sn have been presented and dis- 
cussed in the recent past [|145| - |147| . In ref. [|146|| correlation techniques were 



used in order to determine the multiplicities of H and He isotopes emitted by 
the excited fragments produced at the disassembly stage of an equilibrated 
hot source. The determination of LCP multiplicities and kinetic energies led 
to the fragment excitation energies. These energies per particle are the same 
for all IMFs, which is interpreted as a sign for the fact that thermodynamic 
equilibrium is established at the stage where the source decays. A very de- 
tailed analysis of the same reaction was presented and discussed in ref. ||147| . 
It is shown there that the experimental fragmentation data can be interpreted 
in terms of unique source whose existence could be ascertained by means of 
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kinematical arguments related to global observables like the collective flow 
angle. For events which are selected under specific criteria the angular distri- 
butions of light particles and intermediate mass fragments show the charac- 
teristic behaviour of statistical evaporation, i.e. approximate symmetry with 
respect to 90°, see Fig. || and excitation functions decrease exponentially as 
a function of energy. The experimental results were further confronted with 
simulations by means of the SMM model in its microcanonical version. This 
confrontation confirmed that the assumption of thermodynamic equilibrium 
is compatible with information extracted from experimental data. The anal- 



ysis relies on the so called backtracing method developed in ref. ||110|| which 
solves the inverse problem by tracing back the issue of the process to the 
characteristics of its origin, i.e. the characteristics of an excited equilibrated 
source. This method has also been applied to experiments performed by the 
ALADIN collaboration at the GSI for reactions of l97 Au on C, Al, Cu and 
Pb targets at 600 MeV- A energy p8Hl50f- 



There exist some indications for equilibration from purely experimental 
origin ||151|| , corresponding to so called spectator reactions, i.e. more or 
less peripheral heavy ion collisions, in particular 197 Au on 197 Au at high 
bombarding energies, 0.6 and 1 GeV-A. In this type of reactions it is assumed 
that part of the nucleons, the spectators, form a highly excited system which 
experiences only disordered motion and from which any coherent collective 
motion is practically absent. Several arguments lead to the conclusion that 
these systems could be in thermodynamic equilibrium. Recently the relative 
velocity correlation between LCPs and fragments has been used in order to 
reconstruct the size and excitation energies of the so called primary fragments 
which are produced at the early stage of the collision of Xe on Sn at 32 
and 50MeV ■ A [ |152|| . The constancy of the excitation energy imparted to 
the system with increasing bombarding energy suggests that thermodynamic 
equilibrium may be reached at freeze-out. 

The fragment distributions are invariant with respect to the formation 
process, i.e. the entrance channel of the reaction. As we shall see in section 
5, there has been an attempt to fix a temperature to the spectator system 
by means of the measurement of the isotope ratios of light nuclei (He, Li) 
which are generated through fragmentation |153| | . These ratios are invariant 
with respect to the bombarding energy. The result is however, at first sight, 
not compatible with other observables from which it is in principle possible 
to extract a temperature. This is so, in particular, for the slopes of the 
energy distributions of fragments which, if they are interpreted as Maxwellian 
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distributions, lead to much higher temperatures. Several arguments which 
try to explain this fact have been proposed [151 1 , among them the role played 



by the Fermi motion ||154| of nucleons related to a fast fragmentation process. 
The equilibration problem has also been raised by the FOPI collabora- 
tion in the study of reactions involving more central heavy ion collisions |89 



In this type of events, one observes a collective flow of particles. This fact 
complicates the experimental analysis and does not allow for a clear-cut con- 
clusion about equilibration. Indeed, different models, both time-independent 
and time-dependent seem to be able to reproduce the data ||155|| . Recently 
isospin arguments have been used in order to test equilibration ||156|| . They 
indicate that equilibration is not reached as far as isospin degrees of freedom 
are concerned. 
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4 Percolation models and fragment size dis- 
tributions 

The detection and identification of the charge and the mass of fragments 
generated through violent collisions is of fundamental importance in the 
study of excited nuclei. Many detectors like ALADIN, FOPI, INDRA, the 
EOS setup, the MINIBALL, LASSA, CHIMERA are able to detect and iden- 
tify particles and fragments event by event and to determine their asymptotic 
kinematical properties which constitute the information which can be exper- 
imentally collected. 

As already seen above, the information contained in the outcome of pro- 
cesses like fragmentation are of extreme complexity and the formation mech- 
anism largely unknown, certainly different under different physical conditions 
such as the energy range and the impact parameter. On the other hand, the 
use of the concepts of statistical mechanisms has shown that paradoxically 
the most complex systems can be described in terms of minimum information 
theories which bypass the detailed knowledge of the underlying microscopic 
systems. Here it may be tempting to follow this philosophy and to try to 
describe the outcome of nuclear fragmentation processes within the simplest 
possible physical framework. This has been done starting in the early 80's by 
means of minimum information approaches like percolation models. These 
models ground on purely topological and statistical concepts, they avoid the 
explicit introduction of Hamiltonians. If they work, they raise of course the 
question why and how this is the case. 

In percolation theories, space contains empty regions and parts which are 
occupied by individual objects which can be considered as bound to form 
larger entities defined by means of different criteria. The essential concepts 
are connectivity, localization and percolation which characterizes the prop- 
erties of clusters (animals, fragments). Percolation itself is established when, 
in an infinite space, at least one cluster of infinite size is present. 

Percolation models have been used in many fields which deal with ge- 
ometry, disorder, statistics and physics such as aggregation, localization, 
diffusion, conductivity . . . Before going over to its application in the field 
of nuclear fragmentation it may be worthwhile to recall general concepts, 
definitions and properties which appear in the framework of this type of 
descriptions. 
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4.1 Random-cluster models 
properties 



general concepts and 



The first percolation model was introduced by Broadbent and Hammers- 



ley ||157|| in order to describe the spread of matter through a medium such as 
a liquid through a porous medium. It was then soon recognized that this type 



of model is acquainted to already known models like the Ising model ||158| and 
other statistical spin models like those developed by Ashkin and Teller 
and Potts fll60| . 



All these descriptions can be studied in a general framework called 
random-cluster model. This has been worked out by Fortuin and Kaste- 
leyn ||161 , 162fl . In the first paper |161|1 the authors showed the close connec- 
tion between the connectivity problem of different systems with uncorrelated 
bonds between sites which are occupied by single objects (particles) and the 
well known Ising model. Since Ising systems exhibit phase transitions in 2d 
and 3d spaces it is not surprising that the random-cluster model shows itself 
a phase transition in the infinite space limit. In later developments [11x2 
the authors studied the general topological properties of the aforementioned 
models by means of graph theory, in particular those of the random- cluster 
model. They showed why the models quoted in refs. |158| - |160|] as well as linear 
resistance networks and percolation models can be studied in a unique frame- 
work by means of graphs defined in terms of vertices (sites), edges (bonds), 
connections between edges and vertices, and probabilities p for edges to exist 
or not. In practice they introduced a generalized partition function (cluster 
generation function) and "free energy" which enables to work out physical 
observables like thermodynamic quantities in the case of the Ising, Ashkin- 
Teller and Potts models and correlation functions. 

Further analytical developments on percolation models were pursued by 
Kunz and Souillard ||163| , |164| related to previous work by E.H. Lieb. In the 
case where sites do not interact, i.e. are not correlated, they introduced the 
free energy 
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where p is the bond probability, q = 1 —p, P n the probability that a given site 
belongs to a cluster with exactly n sites, C are all finite clusters containing 
the origin of space 0, \C\ is the size of the clusters and d\C\ the size of 
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the boundary of clusters, i.e. their perimeters (surfaces). The study of the 
properties of f p (h) shows that this function is analytic for h = and p < p c , 
and develops a singularity for p > p c if the system is infinite, indicating a 
phase transition of second order at p — p c . For p > p c the system percolates 
in the sense defined above, i.e. at least one cluster of infinite size appears. 

4.2 Percolation models : definitions and general prop- 
erties 

Before discussing the applications of percolation concepts to nuclear frag- 
mentation we present here some essential properties of percolation models we 
shall need in the sequel. 

Percolation models consist of a discrete number of sites which cover a 
finite or infinite part of a d-dimensional space [|165| - |167| . Site percolation 



corresponds to the case where each site is either occupied by one "particle" 
or is empty. Bond percolation corresponds to the case where each site is 
occupied, neighbouring sites being bound or not to each other. The number 
of neighbouring sites is fixed by the connectivity which depends on the ge- 
ometric structure of space occupation (f.i. a triangular lattice for d = 2, a 
cubic lattice for d — 3, . . . ). The combination of site occupation and bond 
formation leads to hybrid site-bond models. 

Except for the results presented in section 4.1 above [nTT-IIPl], very little 



is known analytically about these models. Most information has been gained 
by means of numerical simulations. 

In the case of site percolation, one fixes a probability p G [0, 1]. For each 
site, one draws a random number 77 taken from a uniform distribution in the 
interval [0, 1]. If 77 < p the site is said to be occupied, if 77 > p it is empty. 
Hopping over all sites, one generates an ensemble of occupied and empty 
sites. If the number N of sites gets very large, one finds Np occupied and 
N(l — p) empty sites. If neighbouring sites are occupied they belong to the 
same cluster. Space is finally covered with clusters and empty space. 

A similar procedure works in the case of bond percolation. For a fixed 
bond probability p consider a pair of neighbouring sites and draw a random 
number rj G [0, 1]. If 77 < p the sites are linked by a bond, if 77 > p they are 
not. Testing this way all possible bonds between neighbouring sites generates 
linked clusters of bound sites, each cluster being made up of an ensemble of 
connected particles. The determination of the whole cluster network will be 
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called event or realization. In each realization there are clusters of different 
sizes. They appear with a given multiplicity which depends of course on 
p. Averaging over many realizations leads to a characteristic cluster size 
distribution in the infinite system. 

Clusters with a fixed size can have different shapes for d > 2. The average 
number of clusters of size r per lattice site taken over a set of realizations is 
given by 



n r 



J29rsP r (l-pY 



where g rs is the number of clusters with size r and perimeter s, the perimeter 
being the number of empty (non connected) sites which surround the cluster. 
Generally g rs is not analytically accessible. We shall come back to this point 
below. 

Individual clusters of size s can also be characterized by a "gyration" 
radius Rn, defined by 



where r*o = 2_, 



i=l 

s 



i=l 



is the centre of mass of the cluster. 

Then the average squared distance between two cluster sites is 



2^ = E 



-> 12 



s 2 



If g(r) is the probability that a site which is a distance r apart from an 
occupied site belongs to the same cluster, one can define a correlation length 



£ such that |167|| 
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2R 2 S is the average distance between two sites of a cluster and sn s the proba- 
bility that a site belongs to a cluster of size s. Since such a site is connected 
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to s other sites the correlation distance £ can be rewritten as 



ta 2E,^. 



e 2 = y f ( 5 ) 

As it has been mentioned above [|161| - |164| there appears a percolation (second 
order) phase transition for d > 2 and some value p = p c < 1- For p > p c 
one or several infinite clusters are present. The system is said to percolate. 
The value of p c depends on the dimensionality of the system, the type of 
percolation (site or bond) and its connectivity (the number of maximum 
occupied sites in the next neighbourhood of a site or possible bonds of a site 
to the nearest neighbours) | |167| ] . 

The central observable in percolation physics is of course the cluster size 
distribution. In practice, nothing is known analytically about it, except 
for conjectures which agree very precisely with numerical simulations. For 
p < p c and clusters of large size (s — > oo) Inn, oc — s, i.e. one observes an 
exponential tail when space occupation (site percolation) or the number of 
bonds ( bond percolation) is sparse. For p > p c and s — > oo \nn s oc — s l ~ l l d . 
As one can see the space dimension enters in the case where space occupation 
or the number of bonds is large. This is not the case in the non percolative 
regime. Finally for p ~ p c and large clusters it has been conjectured [ |167| 
that 

7i, (p) ~ s~ T f({p-p c )s a ) (6) 

which is true for all percolation models, r and a being universal exponents 
and / a model dependent function. For p = p c , /(0) = 1, leading to a power 
law behaviour which intriguinly reminds Fisher's expression (p. As stated 
above, these relations are very accurately verified by numerical tests. 

4.3 Moments of the cluster size distribution and crit- 
ical exponents 

For infinite systems it is possible to work out the expression of the cluster 
size distribution n s (p) in the vicinity of the critical point p = p c . 
By definition the kth moment is 

™<k(p) = ^s k -n s (p) 
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Going over to the continuum limit 

poo 

m k {p) = / dss k n s (p) 
Jo 

and introducing the analytical expression (||) with 

f((p ~ Pc)s a ) = exp(-|p - p c \ 1/a s) 
for p in the vicinity of p ~ p c leads to 

m k (p) ~ Ip-Pcl* 1 " (7) 

t — k — 1 
where //& = (8) 

For k > t — 1 one sees that the moments diverge at p — p c , a sign for the 
existence of a phase transition. Hence, if r < 3 which is the case for d — 2, 3 

m 2 {p) ~ |p-p c r 7 

with 7 = (3-r)/o- (9) 

For d = 3 in a cubic lattice 7 = 1.74. 

Critical exponents characterize other quantities of interest. If in the in- 
finite system P(p) is the probability that a given site belongs to the infinite 
cluster for p > p c , then probability conservation implies 



P(p) + (l-p)+^2sn s {p) = 1 



The second term corresponds to the probability that the site is empty and 
the third that it belongs to a finite cluster of size s. 

At p = p c P{p c ) = 

and Y,s sn s(Pt) = Pc 



Hence in the vicinity of p c , P can be written [ |167|| 



P = YfrM - n,(p)]s + Q(p - p e ) (10) 
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If one introduces the parametrization n s (p) oc s T exp(— cs) with c = at 
p = p c and goes over to the continuum limit with ( |I0|) it comes 



P OC (p - Pc) 13 

where (3 = (r - 2)/<r (11) 

which consistently tends to zero when r > 2 at p = p c . 

Since percolation systems show a continuous transition at p = p c one 
expects that the correlation length £ defined in (Q) shows a singular behaviour 
at the threshold where an infinitely large cluster merges. This is indeed the 
case 

Z<x\p-Pc\-" (12) 

and the positive exponent v can be determined numerically. 

Percolation clusters are generally fractal objects. The radius of gyration 
Rd is related to the volume s through 

s oc Rp 

where df is an effective, not necessarily integer space dimension. From the 
expression (|D it follows that £ 2 behaves like s 2+2 ^ df . The same type of 
derivation as the one used for the moments m},(p) leads to a power law 
expression of £ 2 whose numerator diverges with an exponent (3 — r-\-2/df)/a. 
Hence £ 2 itself diverges like 2/dfa and an identification with (^) leads to 



ov 



1/df (13) 



At p = p c R s oc s 1 '^, hence the size of the largest cluster increases like 
L oc s l l d s . If P is the strength of the infinite cluster as defined above, PL d 
is the volume occupied by the largest cluster. If the linear extension of this 
cluster is of the order of the correlation length £ 

PL d oc L d f 

L oc £ oc \p-p c \~ u 

Recalling the behaviour of P in the neighbourhood of p c one gets 

d f = d-(3/u (14) 
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In summary, general considerations lead to six critical exponents 
r, a, j3, v, df, and 7 which are linked by four relations (|9]), (pTTj) , (|T3| ) and (p!4|). 



Up to now we considered essentially the properties of infinite systems and 
concentrated on their behaviour in the vicinity of the percolation threshold 
which allows to determine critical exponents. These exponents are very use- 
ful quantities since they fix the universality class to which a critical system 
belongs. If excited nuclei are the finite counterparts of critical nuclear matter, 
the singular behaviour of observables at the phase transition is quantitatively 
characterized by critical exponents. We shall try to see and to understand 
below if, how and why percolation concepts could be related to nuclear frag- 
mentation. But nuclei are finite systems. Before we present and discuss the 
models which have been proposed, we present a study of the effects generated 
by the finiteness of percolation systems. 

4.4 Finite size constraints on random-cluster systems 

Percolation deals with the critical behaviour of systems in the limit of 
infinitely extended space. Microscopic and mesoscopic systems like nuclei, 
atoms, molecules, aggregates are finite. As far as percolation concepts can be 
applied to their description and directly compared to experiment as we shall 
do below, it is necessary to get some insight into the behaviour of percolation 
models in finite space. Finite size effects enter because of the presence of a 
surface which induces constraints on the system, especially when the space 
occupation is large. We concentrate here on site percolation. 

The aim of the study concerns the determination of the fragment size 
distribution of a system with V sites, with an occupation number equal to 
or 1 for each site. The number of occupied sites A is fixed, the density is 
p = A/V. Following percolation concepts one qualitatively expects a set of 
small clusters when p is small, hence a cluster size distribution which decays 
monotonously as a function of the cluster size and for large values of p one 
or several large clusters together with small clusters showing a kind of U- 
shape. For intermediate values of p small and intermediate size clusters may 
coexist. These shapes are in fact observed in the experiment. The crucial 
point concerns the behaviour of the distribution for these intermediate values 
of p. 



The behaviour of finite percolation systems has been investigated in |168 
One considers a finite volume V in a d- dimensional space. Space occupation 
by A particles is defined by O = (ii, . . . ,1a) where h, . . . ,%a label A oc- 
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cupied sites. There are ( A ) = V\/A\{V — A)\ possible arrangements. To 
each occupation can be attributed a canonical weight, exp[— /3H(0)] where 
H is the Hamiltonian of the system and (3 can be interpreted as an inverse 
temperature. Define 



H(0) 
with Hij 


= H^ 


A simple choice would be 




H {j = 


- —eLij 



where e is a positive constant and L^ (i ^ j — 1, . . . , V) taking values and 
1 defines a link structure which decides whether two sites are connected or 
not. 

The link structure introduces the concept of clusters as sets c of connected 
occupied sites c = (ii, . . . ,i a ), with size a(c) and a surface (called perimeter 
in section 4.1) made of those unoccupied sites which are the nearest neigh- 
bours of occupied ones. 

Given a cluster of size a and surface s, we call f(c) = V — a — s the free 
complement of c which is such that 

cusu/ = y 

The central observable is the average multiplicity of clusters of a given size 
in the canonical ensemble, allowing for all possible occupations with fixed a, 
weighted by the Boltzmann factor exp(— f3H). 
The average of an observable Q is defined by 

< Q >= \ Yl fi (°) exp[-/3#(0)] 5(a(0) - A) 
ocv 

where obviously 

Z = Z{V, A,f3)=J2 exp[-/3if (O)] 5(a{0) - A) 
ocv 

is the canonical partition function. In particular, the energy U and the 
entropy S read 

U = <H> (15) 

S = \nZ + f3U 
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Notice that, as generally admitted in random-cluster models, there exists 
no cluster-cluster interaction. This is due to the fact that the interaction is 
restricted to nearest neighbour occupied sites. 

The probability to find a cluster c of size a in a configuration is given by 

P(c) = e~ pH ^Z(f, A -a, P)/Z(V, A, (3) 

Then 

<n>=J2 p ( c M c ) 

c 

for any observable Q. 

In particular, the multiplicity of clusters of size a is given by 

m(a) = y^ P(c)S(a — a(s)) 

c 

The essential constraint on m(a) is fixed by the fixed number of particles 

A 

^am{a) = A (16) 

a=l 



Monte Carlo simulations ||168|| reproduce the expected cluster size distribu- 
tions which were predicted above, in particular the £7-shape at large density. 
This is true both in the case when the coupling constant e is zero or finite 
except that the transition from the falling distribution for small p to the 
[/-shape is located at different values of p for different values of e. Hence 
in order to study the consequences of the constraint imposed by (|16|) one 
may restrict this study to the case e = which allows for the use of simple 
analytical arguments. Then 

Z(V,A) = Q 

where s is the surface defined above and g(a, s) the number of clusters of size 
a and surface s. 
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In the case where the system is dilute (p = A/V <C 1) 

J29(a,s)p a (l-pY (17) 



m(a) 



The degeneracy g is a complicated function of the link structure L. There 
exists some knowledge about it [ |165| , |166| , |169|| . In fact g(a, s) differs from zero 
only in a narrow range of surfaces. 
It is convenient to introduce 

g(a,s) =7(a)P a (s) 

where P a (s) is the distribution of surface sizes for fixed a. The surface s 
varies between a minimum and a maximum limit for fixed a. The minimum 
is given for the more compact forms 

s min (a) oc a x ~ ljd for a < — 

oc {V-a) 1 - 1/d for a > — 

The maximum corresponds to chain-like clusters, hence 

s max (a)(x2 + 2(d-l)a (18) 

in a d- dimensional hypercubic lattice. Of course s < V — a. 

The overall behaviour of the surface as a function of the cluster size is 
shown in Fig. |7j. If one defines an average surface 



o\a) 



E sP ^' 



then for the small mass sector one may parametrize a (a) in the following 
form 

a(a) = Aa + o"o aC^o, X < 2d — 2 

and for the large mass sector 

a{a)=K{V-a) a > A , K < 1 (19) 
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where A Q is the solution of s max (A ) = V — A , see flTB|). Hence 

2(d - l)i4o + A = V - 2 

and 

A) = (V-2)/(2d-l) 

Consider first the small mass sector. Since the distribution P a {s) is not 
known, one might investigate different types of distribution. Analytic results 
can be obtained for a Poisson or a sharp distribution. In the simplest case 

P a (s)=5(a(a)-s) (20) 

The distribution may indeed be rather sharp when a does not get too large, 
see Fig. [T[ Then, using ( p0[ ) one obtains 

G(a) = Va(l)[a(a) + a - 2]!/[o! cr(o)!] 

For a, er ^> 1 

G(a) oc a _2 ' 5 exp(a c a) 
with a c = ln[(A + l) A+1 /A A ] 

A similar expression for G(a) is obtained if P a (s) is a Poisson distribution. 
As a consequence 

m(a) oc a~ 2 ' 5 exp{(a c — a)a} 
with a = ln[l/p(l - A) A ] 

a is minimum for a = a c at 

Pc = 1/(1 + A) 

Then m(a) oc a -25 for p = p c . One again finds the manifestation of a power 
law behaviour, which is realized in the case of finite systems. 

If V is large a(a) lies close to 2(d — l)a + 2, hence A = 2d — 2 for cubic 
lattices. Then 

Pc = l/(2d-l)=A /V 
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In fact p c is determined by the crossing of the two topological limits fixed 
by the small and large cluster sector. Hence in finite systems the onset of a 
"critical" regime is correlated with a topological constraint which, in a finite 
but large volume, corresponds to cluster surfaces with a maximum in the 
vicinity of p c V = V/(2d - 1). 

In the case of large p it is sensible to parametrize the surface following 
(|j~9l ) above. Then 

m{a) ocq A - a /{A-a)\ 

with 

q = {(V- A)/ (A + K(V - A))Y- K A(l - K)K- K 

This distribution is monotonically rising with a if K < K cr ~ A/(A+ 1) and 
reaches a maximum at a = A. In all other cases it peaks at q ~ A — a. For 
s m in(a) > V — a, m(a) is of course zero. 

In summary, simple considerations lead to a different behaviour of cluster 
size distributions in different density sectors of finite systems. This behaviour 
is essentially governed by the mass conservation constraint (|l^) imposed by 
the finite volume in which the system is located. For small densities the size 
distribution decreases essentially as an exponential. There exists a density 
p c for which the distribution is a power law or close to it, with a fixed index 
irrespective of the dimensionality of space. For p > p c one observes both a 
decreasing light cluster sector and a large size sector building up. 

The present results are general. As we shall show below they are in quali- 
tative agreement with explicit simulation results and experimentally observed 
cluster size distributions. 

4.5 Models related to percolation and first confronta- 
tions with the experiment 

The first applications of percolation concepts to the description of nuclear 



fragmentation were proposed by Bauer and collaborators [ 170 , 171 |, Campi 
and Desbois ||172 . 



The starting point of Bauer et al. was the observation of the power law 
distribution of intermediate mass fragments flTRpiLKIJ obtained by means of 
fragmentation processes 

A p + A t -»• A f + X 
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where A p , A t , Af are the mass of the projectile, target and a specific fragment 
and X the rest of the reaction outcome. 

In the absence of any precise information about the reaction mechanism 



the authors |170| , |171|| suggested a minimal-information description, following 
the same spirit as the one which prevailed in ref. [|1 73f| . They introduced a 
nucleus lattice model (NLM) in which a set of particles occupy lattice sites 
in a finite space. The presence of bond clusters is tested by means of a 
bond percolation algorithm as described above applied to a 3d cubic lattice 
model in a finite volume. The so called standard percolation theory (SPT) 
deals with infinite systems. The physical framework here is an A + 1 particle 
system obtained by means of energetic collisions of a proton on a nucleus. 
The penetration of the proton generates a cylindrical tube of radius r with 
very excited matter (fireball), in an initial spherical nucleus at an impact 
parameter b. The nucleons of the fireball are kicked out, the remaining 
nucleons are sitting on lattice sites and experience a strong disturbance which 
reflects in the definition of a bond breaking probability p characterizing the 
link between the remaining spectator nucleons. In this description the exact 
number of particles is not fixed and there is no information about the energy 
imparted to the system. The bond breaking probability p is supposed to be 
the larger, the larger the kinetic energy of the incoming projectile proton. 
A Monte Carlo percolation algorithm applied to the system for different 
values of p G [0, 1] leads to different shapes of the fragment (cluster) size 
distribution, ranging from an [/-shape for small p to a typical exponential 
fall-off with no heavy fragments available for large p, as predicted in ref. |168| 



For some p ~ p c ~ 0.74 the width of the mass distribution gets maximum 
and the mean yield shows a power law behaviour for light and intermediate 
mass fragments, with an exponent r ~ 2.21. This value is somewhat lower 
than the one corresponding to the infinite system. 

This bond-breaking model was confronted with experimental results con- 
cerning the reactions p + Ag, p + Ta and p + Au, for proton energies of 
several GeVs. If the bond probability p is parametrized in terms of the 
impact parameter as 

p(b)=p /[l + exp((b-R)/d)] 

where R is the radius of the system, b the impact parameter an d a diffusion 
coefficient, one finds a very good agreement between calculated and exper- 
imental fragment size distributions, see Fig. || by sampling over the whole 



48 



range of impact parameters, b = to R. In order to get a closer link between 
percolation calculations it is sensible to replace p which is not an observable 
by the average over realizations of the fragment multiplicities m obtained for 
a fixed value of p. This allows to relate r extracted from a power law fit 
of the fragment size distribution to an average multiplicity at the "critical" 
point, < m > c . 

The strikingly good agreement between the application of simple per- 
colation concepts and the outcome of inclusive observables extracted from 
complex collision processes raises several questions. The model does not care 
about mass and energy conservation, possible geometrical deformation effects 
when the fragmenting system is generated through heavy ion collisions. The 
bond probability is necessarily related to the microscopic dynamics. It may 
be asked why this unique parameter parametrizes so perfectly the underlying 
interaction between the nucleons in the excited nuclear medium. We shall 
come back to this point in section 5. 



In a similar spirit, Campi and Desbois ||172| developed a continuum site 
percolation approach. A T particles occupy an ensemble of points {r^p^i = 
1, . . . , At} in phase space chosen in a classical one-body Wigner distribu- 
tion fwifiP) determined by a one-body Hamiltonian with a Saxon- Woods 
potential. The collision process leaves pA? nucleons of the original nucleus, 
(1 — p)A T nuleons are kicked off the original nucleus, p G [0, 1] depends on the 
characteristics of the reaction like impact parameter, energy, target size . . . 
Clusters are generated by means of phase space prescriptions which decide 
whether particles belong to fragments or not. A given number of particles 
forms a bound cluster if 

dij = \n -fj\- \pi -Pj\ < 2.5H 

for every pair of nearest neighbours belonging to the supposed cluster, where 
the maximum action on the r.h.s. of (|4.5|) corresponds to the value obtained 
for the deuteron. This condition is supplemented by a compactness constraint 
which imposes that the mean square radius Ra and mean square momentum 
Pa verifies 

Ra < (H-e)roA 1 / 3 (21) 

where p F is the Fermi momentum and e ~ 0.12. If the constraints ( pT|) 
are not fulfilled the longest link d^ is cut and the determination of clusters 
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PA < (1 + e) [■=) P F 



starts again. The first condition takes care of the balance between volume 
and surface energy of the cluster, the second eliminates nucleons with too 
large momenta. The model is constructed with percolation concepts which 
do however not correspond to the standard ones, essentially because the 
generation algorithm is not fixed by a single and simple bond probability 
parameter. The cluster multiplicity obtained in this way can be fitted by 
means of parametric analytical expressions [|172|| . The generated clusters 
are excited objects whose excitation energy E* can be estimated. Particle 
emission from fragments and their fission are taken into account. 

In order to enable the confrontation of the model with experiment, the 
quantity p which fixes the number of remaining nucleons after the collisions 
is related to the number of nucleon-nucleon collisions which are experienced 
in the system for a fixed impact parameter b and nucleon mean free path 
A which varies with the energy. This leads to a distribution for p, p(p), 
for p larger than a value p m in corresponding to b = 0. Depending on the 
physical conditions, p can reach or not the percolation threshold. One can 
then calculate explicitly inclusive mass yields for fragments of size A F 



a(A F ) = dp p(p) m(p, A F ) 



where m(p, Ap) is the multiplicity of fragments of size A F , for a fixed value 
of p, and the total reaction cross section 

** Pmin 

The model reproduces very well the fragment size yields u{A F ) for low and 
high energy data corresponding to 20 Ne + 181 Ta at 400 MeV-A, p + 181 Ta at 
5.7 GeV and p + 181 Ta at 340 MeV. In the lower energy data one observes 
the opening of the fission channel which is predominant along with particle 
evaporation at this energy. For high energy reactions, p+Xe at 80 — 350 GeV, 
both the slope and absolute value of a(A F ) are well reproduced. The trend of 
the apparent "critical" exponent a(A F ) oc A F T where E is the bombarding 
energy shows a minimum r ~ 2.4 at E = 50 GeV. This is understood in the 
following way. When E is small p m in is large, p c cannot be reached, hence 
the system is essentially made up of one large cluster and a set of small 
particles. When E increases p min decreases, p c can be reached and one gets 
a power law distribution over the whole range of fragments sizes, except for 
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finite size effects. For large value of E p min continues to decrease, more 
and more smaller fragments are generated and the apparent slope of the 
distribution increases again. At very high energy the nucleon mean free path 
X(E) becomes constant, hence p m in and r saturate. 

4.6 Relevant observables and percolation analysis of 
experimental data 

The first extensive confrontation between standard percolation theory 



and nuclear fragmentation was initiated by Campi [|174|| . The phenomenon 
was interpreted as a bond-breaking process between bound nucleons due to 
the energy supplied to the system under violent collision conditions. As 
discussed in section 4.3, the moments m k of the fragment size distribution 
get singular at p = p c for k > 2, so do the normalized moments S k = m k /mi. 
Hence one expects that they show a maximum in the vicinity of p ~ p c . Since 
m k oc exp(— fi k ) with \x k = (r — 1 — k)/a (see (|7|) and (|^)) at p c , there should 
exist a linear relation between In S k and In Si at p c . For k = 2 and I = 3 the 
slope can be written as hi S3 = AlnS^ with A = 1 + l/cry and one can in 
principle read off the value of the power law exponent r. 

The introduction of m k and S k allows for a direct comparison between 
the moments obtained from bond percolation models (here 3d cubic) and 
those obtained from the experiment, m k , which can be determined from the 
fragment size distributions of a set of events labelled by the index j by taking 
an average over all available events. Another interesting observable is given 

by mg 



■n = ^ (22) 

which is related to the variance a of the fragment size distribution by 

1 i 2 2/2 

72 = 1 + a ■ ■m Q /m 1 

and should also show a sizable enhancement at the "critical" threshold. 
Finally the last observable proposed in ref. ||174|| concerns the behaviour 



of the largest cluster A max in the system. Indeed, in the infinite system, 
A m ax should get infinite at criticality, hence it is expected to get large when 
5*2 gets large. Hence the correlation between A max and S 2 is expected to 
contain information of interest. 
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Bond percolation calculations in a cubic lattice in a finite system were 
performed for different sizes and S 2 ,S 3 , 72 and A max were calculated. The av- 
erage quantities were confronted with the outcome of a few hundred events 
generated by the fragmentation of l97 Au nuclei in emulsions ||1 76| . Typi- 



cal fission events involving two large fragments were excluded. Results are 
shown in Fig. [^ which relates 5*3 and 5*2 for both experimental events and 
percolation simulations. "Critical" events correspond to large values of S2 
and S3. Violent (p = 0) and smooth (p — 1) collision events correspond to 
small values of S2 and 53. One observes that there is a striking agreement 
between the two results. It should be noticed that in both cases one mixes 
events which correspond to different values of p (as far as calculations are 
concerned) and multiplicities which are related to the degree of violence of 
the process (as far as the experiment is concerned). They both show the 
expected linear behaviour which can be related to the presence of "critical" 
events. It is also remarkable to notice that the slopes A ~ 2.22 are in both 
cases very close to the value A = 2.25 corresponding to bond percolation 
in the infinite system. The same type of agreement can be observed when 
dealing with S 5 vs. S2 ||1 74|| . In Fig. |10| 72 is drawn as a function of the 



average multiplicity tjiq. The similarity of simulation results and experiment 
is again remarkable. One observes an enhancement in the neighbourhood of 
m — 0.25 which, when related to the bond probability p, corresponds to 
p ~ p c in 3d cubic systems. The fact that 72 > 2 for m ~ 0.25 indicates 
a power law distribution. Finally Fig. [TI| compares the behaviour of A max 
as a function of S2 in both cases. "Critical" events are concentrated in the 
right corner of the figure. One observes again a strong agreement between 
calculations and experiment. 

The possibility to extract values of A and r and the fact that these are 
close to standard percolation results may induce the temptation to conclude 
about the universality class of the apparent critical phenomenon unravelled 
by the present studies. One must however be careful, because different types 
of transitions may lead to exponents which lie close to each other but corre- 
spond to different classes. The determination of exponents needs some care. 
We shall come back to the problem concerning the determination of critical 
exponents in section 5. 
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4.7 Extensive tests on peripheral collision data : anal- 
ysis of ALADIN experiments 

Many experimental efforts have been made over the last decade by the 
ALADIN collaboration of the GSI in order to collect a large amount of sig- 
nificant information about particle and fragment multiplicities generated in 
the fragmentation of the projectile by means of energetic peripheral heavy 
ion collisions. Such precise and complete information allowed for a care- 
ful comparison with theoretical approaches, in particular percolation and 
percolation-inspired models. The first measurements which were performed 
established correlations between the multiplicity of light particles Mi p , the 
charge Z max of the largest fragment in an event and the sum of the charges of 



the fragments with Z > 2, Zf, und ||148|. It was found that the correlation be- 



tween the multiplicity of intermediate mass fragments (IMFs), i.e. fragments 
with 3 < Z < 30 and Z bound was independent of the target nucleus. This was 
interpreted as a sign for the equilibration of the projectile remnant (spectator 
particles) before decay. The first detailed investigation of the average multi- 
plicity of IMFs, < Mimf >, showed that with increasing energy transferred 
to the projectile in reactions of Au on C, Al and Cu targets < M IMF > first 
increases to a maximum and then decays [ |149|| . This could be interpreted 
in the following way. For small energies peripheral collisions lead essentially 
to one large projectile remnant and a few light particles. With increasing 
energy, the projectile breaks up into more and more IMFs. If the energy gets 
very large most fragments are very small, (Z < 2) and hence < Mjmf > de- 
creases again. These results were then compared to theoretical investigations 
of < Mimf > vs. Zb oun d ||150| by means of the statistical models of Copen- 



hagen [|123| - |125|1 and Berlin [|122j| , and a sequential decay scenario ||177|| . It is 



seen that sequential evaporation cannot reproduce the data. 

An extensive analysis concerning the reactions with Au on C, Al, Cu and 



Pb at 600MeV-A was performed in ref. |178|| . The charge distribution was 
fitted by a power law and the exponent r showed a minimum somewhat 
below 2 when represented as a function of Zbound- The average largest frag- 
ment charge < Z max > increases with increasing Zbound- Other correlation 
functions like 72 (see above) and 

A 3 = \{Z max - <Z>f + (Z 2 - <Z>) 2 + (Z 3 - < Z >) 2 ] 1/ Vv / 6 
where < Z >= -(Z max + Z 2 + Z 3 ) 
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and Z 2 and Z 3 are the second and third largest fragment charges were deter- 
mined. 

All results showed the striking feature that it is a simple-minded 3d site- 
bond percolation model with an occupation probability p s for the site and a 
fixed bond probability pb = 0.45 obtained through a fit to the Z max distribu- 
tion which reproduces the best the available data. Masses in the percolation 
calculations were converted into charges by means of a semi-empirical mass 



formula [|178|| . The investigation of projectile fragmentation was finally ex- 
tended to other projectiles like 129 Xe and 238 JJ with Be, C, Al, Cu, In, Au and 



U targets and energies ranging from 400 to 1000 MeV -A [|179||. The striking 



feature about the results concerns the energy independence of < M IMF > as 
a function of Z bound and its scaling properties, i.e. all curves get superposed 
on each other if Zbound is rescaled to Zbound/ Z p , where Z p is the charge of the 
projectile as it can be seen in Fig. [T^. Statistical multifragmentation models 
reproduce qualitatively if not quantitatively the experimental results. 

Recent ALADIN data concerning fragment size multiplicities for 197 Au 
on different target reactions at 600 MeV ■ A were analysed by means of a 
standard 3d cubic bond percolation model ||180|| . Depending on the value 
of the impact parameter, the number of spectator nucleons which remains 
in the projectile remnant is of course different. In order to decide about 
this number, a simulation by means of a BUU transport equation of the 
reaction process was performed. The spectator nucleons were arranged on 
a cubic lattice in a compact spherical arrangement and a bond percolation 
algorithm was introduced. Masses were converted into charges by means of 
the relation 

Z = A/(1.98 + 0.0155A 2/3 ) 

for A > 4. A filter was applied to the results of the calculations in order 
to allow for a direct confrontation between the model calculations and the 
experimental data. A whole set of observables was analysed in this way. 
Fig. [T3| shows the average value of A 12 = (Z max - Z 2 )/(Z max + Z 2 ), A 23 = 
[Z 2 — Z 3 )/(Z 2 + Z 3 ) and A 3 defined above as a function of Zbound- As one 
can see the agreement with the experiment is nearly perfect and this is also 
the case for all other observables. 

A new analysis of the data [|181|| confirms this statement. The authors 
also use a standard 3d bond percolation algorithm. The number of nucleons 
Zps which is taken into account ||182|| is parametrized in terms of Zbound- 
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Simulations are done for fixed Zps and observables are determined for fixed 
Z bound- The analysis is again parameter free. The calculated quantities are 
the fragment charge distributions, the average size and fluctuations of the 
largest fragment, and the largest fragment distribution. The confrontation 
also shows that there appears a "critical" regime, characterized by a power 
law charge distribution with r ~ 2.2. In the neighbourhood of the "critical" 
point the distribution of fragments with charge z can be parametrized as 

n(z) = n c (z)f(z/z) 

where z = m 3 /m2, m^ being the moment of order k defined in section 4.3. 

The authors consider two interpretations of these results. The first relies 
on the fact that a fragmenting nucleus is made of particles which interact 
essentially as nearest neighbours through a short range two-body potential. 
If the relative kinetic energy between pairs of neighbouring particles is larger 
than the potential energy they experience a bond is broken and hence in- 
dependent fragments are generated if the excitation energy is large enough. 
The second interpretation requires the assumption of thermodynamic equilib- 
rium. At different fixed densities it is the temperature T of the system which 
fixes its fragment content. For a given density p of the system, there exists 
a temperature T(p) where the distribution looks like a power law. We shall 
come back to explicit classical microscopic model illustrations of both inter- 
pretations in section 5 below. It is clear that thermodynamic equilibrium 
is not a necessary prerequisite for an accurate description of the fragment 
content of energetic nuclear collisions. 

4.8 Comparison with other fragmentation models 

The intriguing and impressive success of classical percolation models in 
the description of finite nuclear fragmentation events and the existence of a 
continuous phase transition in infinite percolation systems triggered many at- 
tempts to compare and possibly relate different models with the percolation 
approach. Such attempts started with the Statistical Multifragmentation 



Model SMM [125]. The authors tried to establish a link between the percola- 



tion bond probability p and the excitation energy available in the fragmented 



system described by SMM ||124| at freeze-out. They fixed p = Vo/(Vo + V/) 
where Vq is the volume of the initial compact system and Vf the volume 
at freeze-out. The investigations revealed that such a link is not clearly es- 
tablished and led to the conclusion that inclusive fragment size distributions 
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may not contain enough information in order to distinguish the purely statis- 
tical features of the fragmentation process from those which concern physical 
information carried in phase space models. 

A detailed comparison of the fragment content obtained in the frame- 
work of the MMMC model |7JJ with percolation was carried out in ref. [ [183 



by means of an analysis of the relation between the charge of the heaviest 



fragment and the second moment of the charge distribution ||174| . Close sim- 
ilarities between MMMC and percolation models were found, in particular 
the "critical" zone which was discussed in section 4.6. The main qualitative 
difference between the two approaches is the manifestation of the presence 
of fission events (two heavy fragments) which are absent in the percolation 
simulations. This was interpreted as being due to the presence of the long 
range Coulomb interaction which is not taken into account by the percolation 
model. Three-body correlations between the three heaviest fragments were 
also in agreement with the percolation predictions. So called cracking events 
with more than two large fragments emphasized the role and importance of 
the Coulomb interaction. 



Moment correlations 173, 1741 have been used in order to test different 



models and possibly discriminate between them [|184]|. The models which 



are introduced are an equal probability model for the generation of differ- 
ent partitions |186| , |187| ], time-dependent sequential decay processes pi] , |185 



and a standard bond percolation model on a finite cubic lattice. All models 
agree qualitatively with percolation calculations but do not reach quanti- 
tative agreement. The peak corresponding to "criticality" in percolation is 
always present, but it is shifted in the multiplicity and its height is not re- 
produced. As expected, the agreement between percolation calculations and 
experiment is excellent. 

A model analysis of data has also been performed by Kreutz et al. [|178 



with a statistical decay model |177|| , SMM and, as already discussed above, a 
site-bond percolation model. It appears again that all models lead to quan- 
titatively different results. Sequential decay fails to reproduce charge cor- 
relation measurements, < 72 >, charge distributions and IMF multiplicities 
come out more or less closely to experiment in SMM, although the average 
largest fragment charge is too small. 

4.9 Final remarks 

If one restricts the information concerning the fragmentation of highly 

56 



excited nuclei to the description of the asymptotic fragment size distribution 
it appears clearly that parameter free standard percolation models are the 
only ones which are able to reproduce the experimental data quantitatively. 
It raises the up to now unanswered question why this is so. It shows that 
a minimum information model is sufficient in order to explain the a priori 
complex nuclear fragmentation process. The agreement is by essence inde- 
pendent of any detail concerning explicit two-body potentials. As it will be 
seen later, simple models governed by short range interactions are able to 
reproduce percolation results. The long range Coulomb interaction which 
induces fission-like fragmentation is out of reach. The fact that its absence 
does not affect the physics of the process may be understood if the excita- 
tion energy is large enough. The influence of the Coulomb interaction will be 
more extensively discussed in the framework of lattice models in section 5. 
Taking a pessimistic attitude, it is tempting to conclude that the success of 
percolation concepts is the sign that the information content of the process 
is poor. But there remains the fact that the result is not trivial since the 
model is not. 

The success of a percolation description raises the problem of the pres- 
ence of a phase transition in the infinite nuclear system, its origin and its 



nature. The formal similarity with Fisher's model |19[ is striking, since even 
the order of magnitude of the critical power law exponent r is reproduced. 
Fisher's model leads to a thermodynamic phase transition. Does there ex- 
ist a connection of the second order percolation phase transition and the 
thermodynamic properties of the system ? 



The outcome of the analysis of finite systems |168|1 should be kept in 



mind. Mass distributions of finite systems are able to show a behaviour 
which corresponds or is close to a power law behaviour, which may simply 
be induced by particle number conservation in a system confined in a finite 
volume, suggesting that a power law behaviour may not necessarily be the 
sign for the existence of a critical phenomenon. 

Experimentally, charge distributions are only a small part of the avail- 
able information on nuclear fragmentation. This information has to be re- 
produced. It has to be related to the properties of observables related to 
the energy, in particular to the thermodynamic properties of the system if 
it is in equilibrium. This point has already been considered in section 4.8. 
It will be raised again in section 5. There we introduce generic microscopic 
cellular and lattice models which are close in spirit to percolation models but 
also hopefully realistic enough to be able to reproduce the thermodynamic 

57 



properties of excited fragmenting nuclei. 
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5 Lattice and cellular model approaches to 
nuclear fragmentation 

Standard percolation provides a surprisingly good description of the 
asymptotic fragment multiplicity distributions generated in energetic nuclear 
collisions. By essence, percolation models cannot describe the thermody- 
namic properties of the decaying systems since they do not explicitly inte- 
grate the properties related to thermodynamic quantities like energy, volume, 
pressure. In the preceding section we presented models which try to link the 
purely geometrical aspects of percolation models to considerations related 
to phase space ||125| , |170H173| , |183Jl . This link introduces phenomenologically 



reasonable but arbitrary parametrizations of the bond probability, generally 
in terms of excitation energy. The procedure may not be totally satisfactory. 
A hint to a more consistent approach can be found if one remembers the inti- 
mate relation between standard percolation models and spin models like the 
Ising model and its generalizations |162| , |163| . These last models introduce 



explicit Hamiltonians and are generic for many physical systems. The afore- 
mentioned formal relation induced the introduction of classical microscopic 
models, starting with Ising itself. These models allow for a thermodynamic 
description and the study of the fragment content of the system in a unified 
framework. They can be considered as a simplified though hopefully realistic 
description of a nucleonic system. 

In the sequel we aim to introduce and discuss different types of cellular 
and lattice models which have been proposed and studied by different groups 
in the recent past. In order to put this study in a comprehensive perspective 
we first recall some fundamental and useful concepts concerning statistical 
ensembles and thermodynamic phase transitions in finite and infinite systems 
at equilibrium. 

5.1 Ensembles, thermodynamic stability, phase tran- 
sitions in finite and infinite systems 

We summarize here some general definitions and concepts which will be 
used for later developments and discussions presented below. They can be 
found in textbooks, see f.i. refs. |188H190| • 



If one considers a large classical microscopic system which is in ther- 
modynamic equilibrium, each particle is characterized by its phase space 
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coordinates {fi,Pi,i = 1, . . . ,N}. If {ri,Pi} is considered as a point in this 
space, there can exist a large number of points which correspond to the same 
macroscopic state. Each set {fi,pi} defines a realization of the system. A 
collection of realizations which correspond to the same macroscopic ther- 
modynamic state is called an ensemble with an arbitrary fixed number of 
physical properties such as a fixed number of particles, energy, temperature, 
angular momentum and possibly other conserved quantities. 

In statistical mechanics there are essentially three ensembles which are 
used in practical applications. The most common one is the canonical ensem- 
ble which corresponds to a closed isothermal system with a fixed number of 
N particles and which exchanges energy with an external ensemble generally 
called reservoir. If the system is made of N t subsystems of rij subsystems 
with energy Ei (i — 1, . . . , Nt) 



J2 n ' 

i 

^2 niR ' 



N t 
E 



(23) 
(24) 



The thermodynamic probability which characterizes a macroscopic state is 

W({m}) = N t \/Y[' l 



n;\ 



and the entropy reads 



k\nW n 



where k is the Boltzmann constant and W max is the maximum of W consistent 
with the constraints (^) and fl24). The average entropy S = Tt/Nt can be 
written as 
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Pi is the weight of the subsystem i and /3 can be identified with the inverse 
thermodynamic temperature T _1 . The knowledge of the entropy allows for 
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the determination of all thermodynamic functions and quantities such as the 
free energy of a system in a volume V 

F(N,V,T) = U-TS= -kTlnZ 

where U is the energy. 

The concept of temperature has been extensively used in the canonical 
framework in low energy nuclear physics, in particular in Hartree Fock models 



23] and compound nucleus theory [73]. Its use implicitly implies that the 



energy is not a rigorously fixed quantity. It raises also the question of the 
concept of temperature in very small systems like nuclei |191|| . 



In the grand canonical ensemble both energy and number of particles are 
only fixed in the average by supplementing eqs. (E3) and (53) with 



j2 n i N i 



N 



where iVj is the number of particles in subsystem i. Then 

S = -A^Tihiy, 

i 

where % = e -M+Pi»«/E, 



£ 



e 



and \x is the chemical potential which fixes the average number of particles 
of the system. One can again work out all thermodynamic functions, in 
particular the Gibbs free energy G = F + PV where P is the pressure in the 
system. 

If the system is closed, with fixed energy and fixed number of particles, 
it is described in the framework of the microcanonical ensemble. The tem- 
perature T is then no longer a natural concept. Mathematically speaking 
the partition function with fixed temperature is the Laplace transform of the 
partition function with fixed energy, hence energy and temperature are conju- 
gate through a Laplace transformation. The temperature can be introduced 
through the thermodynamic relation 

T- 1 = ( £) 

\dEj v 
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where the entropy 

S= -k^PilnPi 

i 

and Pi is the weight of subsystem i. The principle of equiprobability of 
occupation stipulates that 

Pi = 1/fi for E-AE<Ei<E 
= otherwise 

where Q is the partition sum of the system, the number of systems in the 
ensemble. 
Then 

S = k\n£l 

This shows that the central quantity to be known is the partition function 
which counts the number of states at energy E or the directly related entropy 

S. 

5.1.1 Stability of thermodynamic systems 

Microscopic many-body systems can be microscopically unstable. The 
Le Chatelier principle stipulates that any spontaneous change in the param- 
eters of a system which is in stable equilibrium will give rise to a process 
which tends to restore the system to equilibrium |188|| . This leads to the 
inequalities 

C = ^>0 (25) 

K = -^>0 (26) 

for the heat capacity C and the compressibility K which characterize the 
change in the quantity of heat Q when the temperature changes and the 
change of the volume V of the system if the pressure P varies. Eq. (|25|) con- 
cerns thermal equilibrium and (^6|) mechanical equilibrium. These inequal- 
ities show that the free energy F is a concave function of the temperature 
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T and a convex function of the volume whereas the Gibbs function G is a 
concave function of both T and V. Indeed, in particular 

d 2 F\ fdS\ 1 n 

-( — I = ~-C V < 



dT 2 J v \dTj v T 

d 2 G\ fdS\ 1 



dT 2 J p \dTj p T Cp<0 

for fixed volume V and pressure P respectively. 

If (^) and (^) are not realized the considered system is unstable and 
cannot be described in the framework of equilibrium thermodynamics. 

5.1.2 Some essential reminders concerning phase transitions 

Phase transitions can occur in many types of systems. In the traditional 
acceptance of the word they are observed in infinite systems. Only there can 
the correlation length between different points in the system or thermody- 
namic quantities like the heat capacity get infinitely large. In practice the 
phenomenon can be observed and mathematically extrapolated to infinite 
systems in macroscopic large systems like those which are at hand in con- 
densed matter physics where the phenomenon has been studied for all kinds 
of materials. 

There exist different types of phase transitions in nature. They have been 
classified in universality classes. Thermodynamic transitions are said to be 
of first order if there exists a discontinuity in one or more first derivatives 
of the appropriate thermodynamic potential. If, for instance, the entropy is 
not continuous as a function of temperature then the transition is associated 
with latent heat generation. 

If the first derivatives are continuous but second or higher ones are dis- 
continuous the transition is said to be of second order or continuous. There 
appears then divergences in quantities like the susceptibility, correlations get 
infinite and decay as power laws. The order of a transition is a universal 
characteristic of the phenomenon. Different though closely related systems 
can experience different types of transitions. Such is for instance the case of 
the Ising model which shows a second order transition at some fixed value 
T c of the temperature and a first order transition for T < T c for B — if 
the system experiences a magnetic field B. Then the magnetization which 
is the derivative of the free energy with respect to B shows a jump. It is 
continuous at T = T c but shows an infinite slope there. 
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Phase transitions can be characterized by order parameters. These are 
useful physical quantities which are a help in order to distinguish between 
different phases. They can be chosen arbitrarly among the extensive vari- 
ables. In spin systems the magnetization works as an order parameter. For 
each phase there exists a characteristic value of the order parameter, the 
magnetization is zero above T c and finite for T < T c . Order parameters be- 
have either continuously or jump at a transition crossing point. This is the 
reason for their usefulness when trying to find the order of a transition. 

5.1.3 Critical points, singularities, universality and critical expo- 
nents 

Phase transitions occur at specific values or intervals of values of the ther- 
modynamic variables such as the temperature. As already seen in section 4 
these values which correspond to so called critical points, lines, surfaces, 
are marked by divergences in physical observables. In the case of a ther- 
modynamic transition one expects that in the vicinity of a critical point T c 
observables may behave like 

limF(t) = A\t\ x (l+pt Xl + •■■) 

t-i-0 

with Ai > 0, as a function of the dimensionless parameter t = (T — T c )/T c . 

For instance for a fluid system one observes that the specific heat at con- 
stant volume Cy behaves like Cy ~ |t| _a , the liquid-vapor density difference 
{pi — Pg) ~ ( — i) i the isothermal compressibility Kt ~ |t| -7 and the corre- 
lation length £ ~ \t\~ u , where all the indices are positive. The interest in the 
knowledge of critical exponents lies in the fact that they possess a univer- 
sal character. While the location of the critical exponent itself may depend 
on the details of the interparticle interaction, the exponents are universal 
in the sense that they depend only on quantities like the dimensionality d 
of space and the symmetry properties of the order parameter. Evidence for 
this has been verified number of times, in particular for the first time by 
Guggenheim | |192| in the phase separation line of the liquid-gas transition 



for different species, see Fig. [L4L for which it can be seen that the exponent 
j3 = 1/3. The very same value for (3 comes in the 3-dimensional Ising model, 
the magnetization M at T = T c behaves like ~ (— i) 1 / 3 . These properties 
lead to define universality classes for systems showing the same behaviour at 
phase transitions. They allow to classify a priori very different physical sys- 
tems by means of their critical exponents which characterize the behaviour of 
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their order parameters at a critical point. Scaling relations which govern the 
behaviour of physical quantities when one looks at systems at different scales 
allow to establish inequality or equality relations between the exponents. 
This has already been shown in section 4 in the framework of percolation, 
see also below. 

The importance of critical exponents raises the question of their explicit 
determination. Generally this can only be done numerically. Since critical 
exponents only make sense in infinite systems it is necessary to find methods 
which lead to an extrapolation of the behaviour of the critical observables 
to the infinite limit in order to determine them. There exists different, more 
or less sophisticated methods and numerical techniques which allow for very 
precise determinations of exponents, see f.i. refs. 



A whole class of methods relies on the so called finite size scaling assump- 
tion (FSS). The scaling hypothesis implies that in the neighbourhood of a 
critical point a system shows self-similarity properties when it is examined at 
different scales. These properties are reflected in the mathematical behaviour 
of thermodynamic functions and physical observables which show so called 
homogeneity properties, i.e. 

f(ab) = A(a)f 2 (b) 

where fi{a) = a p , p is a real number, and for an arbitrary number of variables 

f(Xx 1 , Xx 2 , • • • Xx n ) = A p /Oi, x 2 , ■ ■ ■ x n ) 
If, for the case of two variables, one chooses A = x^ 1 then 

f{x 1 ,x 2 ) = z?/(l,s 2 /si) = x\g{x 2 /xi) 

which shows that in practice F does not depend on the two independent 
variables x\ and x 2 , but on X\ and the ratio x 2 jx\ whatever x 2 . For fixed xx, 
g is a scale independent function oi x 2 /xi, g{Xx 2 /Xxi) = g(x 2 /xx). The FSS 
assumption allows to relate the behaviour of an observable at some scale to 
its behaviour at another scale. 

As an example of application consider an observable O^(At) which de- 
pends on a variable At = (t c — t)/t c where t — t c is a critical point of 
the infinite system and describes a finite system of linear dimension L. If 
one defines x(L,t) = C,i(t)/L where £l(£) is the correlation length in the 
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finite system for a fixed value of t in the neighbourhood of t c , FSS stipulates 



that |19^ 



L (At) = O(At)Q (x(L,t)) (27) 

where O(At) corresponds to the value of the observable in the corresponding 
infinite system and Qq(x(L, t)) is a scaling function which is universal in the 
sense that it does not depend on t but only on the ratio x in the neighbour- 
hood of t c , i.e. all points corresponding to different t's lie on the same curve 
when Qo is plotted as a function of x. If Ox, (At) converges to a finite O(At) 
with increasing L, Qo(x) can be determined from ( P7|) and close to t c one 
expects that 

O(At) = O L (At)/Q (x) ~ (t c - t)" a (28) 

if O shows a critical behaviour, and a is a critical exponent. If Qq(x) can be 
obtained numerically it is generally fitted to a polynomial in x 

Qo(x) = l + b 1 x + b 2 x 2 + •■■ 

with Qo(0) = 1 if ^i(t)/L goes to zero when L goes to infinity. This type 
of methods has been applied in different places, in particular to percolation 
(t = p) and spin systems [|195 - |197| . Their use leads generally to reliable 



and accurate values of t c and critical exponents. Their practical implemen- 
tation shows however that the numerical determination of these quantities 
can be a cumbersome task. Other methods have been proposed by different 



authors ||55| , |198|1 , in particular in order to fix the exponents which describe ex- 
perimental fragment size distributions in the vicinity of the expected critical 
point where the distribution should follow a power law in the corresponding 
infinite system. However these methods rely on expressions which are valid 
for large systems and do not extrapolate to the thermodynamic limit. Even 
though the information on criticality extracted this way looks reasonable it 
remains suspicious to use it in order to conclude about the determination of 
the universality class of the phase transition which is the ultimate goal of 
such an analysis. On the other hand, it is difficult to believe that the use of 
scaling assumptions implemented by means of FSS can be used in practice 
in the framework of the analysis of experimental events since, among other 
severe problems, the number of points necessary to construct scaling func- 
tions would be rather restricted because of the small interval of sizes which 
can be reached with nuclei. 
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Critical exponents are rather refined quantities. There remains hope to 
be able to detect the existence and eventually the order of phase transitions 
in nuclear matter by other means. Considering thermodynamic aspects of 
fragmentation, this can possibly be reached restricting ones considerations 
to finite systems. The question concerning the order of a phase transition 
has been raised recently on theoretical grounds ||199| , p00(] in the framework 
of the Lee and Yang theory [|201|| and the approach of Grossmann and col- 



laborators |[202|1 . The criterion which allows to distinguish between different 
orders concerns the location of the zeroes of the grand canonical and the 
canonical partition function in the complex temperature plane. But these 
considerations are theoretical, possibly difficult to implement in practice for 
interacting systems and the connection with experiment is lacking. In section 
6 we shall show that there may exist simple operational means to overcome 
the problem. 

Finite size effects can hide or even delude the truth with respect to the 
order of a phase transition. Different ensembles may lead to a different 
behaviour of observables when applied to finite systems. We shall indeed 
see below that characteristic features of observables when calculated in the 
framework of a specific ensemble may be of some help in order to decide about 
the nature of the transition. The assertion has however again to be taken 
with a grain of salt since there does not exist any guarantee that these char- 
acteristic features will maintain their validity in the thermodynamic limit. 

5.2 Lattice models 

We present and discuss here different lattice models which have been 
introduced and studied by different groups. The essential point of inter- 
est concerns the existence of signs for the existence of one or several phase 
transition(s) in nuclear matter. 

5.2.1 The Ising model as a random cluster model 

The success of percolation raised the question for the reasons of its 
achievements. As we saw in section 4, classical lattice systems described by 
the Ising, Ashkin- Teller and Potts models are generally close to percolation 
models |161| , |162|| . They describe Hamiltonian systems, hence they introduce 



an explicit, physical description of an iV-body system. This was the main 
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motivation for the introduction of the Ising model in order to interpret the 
cluster size distributions of an excited, fragmented system |169| , [203|| . 



In order to be able to get a clearer understanding of the behaviour of 
observables the model was first studied in d = 2 which allows the use of 
analytical techniques ||169| . The Hamiltonian reads 



N M N M 



H — — J\ 2_^ 2^ a nm cr n+lm ~ ^2 2_^ /_^ a nm^r, 



n=l m=l n=l m=l 



where J\ and J2 are positive constants and a nm is a classical spin (= ±1) 
located at the position (n, m) on a 2d lattice with N lines and M columns. 
An ensemble of sites are said to form a cluster if the corresponding spins 
show the same value and each site of the cluster has at least one neighbour 
site with a spin of the same value (connexity). The perimeter of a cluster is 
the ensemble of nearest neighbour sites whose spin takes the value which is 
opposite to the spins in the cluster. With this definition and after a mapping 
of the classical system into a quantum system of independent fermions ||204|| it 
is convenient to introduce projection operators along lines (I) and columns(c) 



x ■> 

nm u nm+L 



Pi C) = l(l + rrL,< 



and 



2 



1 



nm n+lm) 



Pi -(1 -a* a, 



X 



2 

,(c) l 



nm nm+lj 



r - — q^ u nm u n+lm) 

where a x is the standard Pauli matrix. The correlation operators P + and 
anticorrelation operators P_ have a classical interpretation, P\_ ,c = 1 and 
P_ = for aligned spins, and P|_ = and Py c > = 1 for anti-aligned 
spins. 

These operators allow to define the probability of presence of a cluster of 
a given size and geometry. 

r(i) r(0 

n« =< p+ + • p- > (29) 
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where the index i characterizes the linked cluster which is considered, L+ 
is the number of bonds between cluster sites and L_ the number of bonds 
between the surface of the cluster and the perimeter sites. For a given cluster 
size Ns, one finds 

2L ( i ) + L W = 4iV s 



The brackets in fl29|) indicate a trace over the many-body quantum states 



im 



The multiplicity corresponding to clusters of a fixed number of sites Ns 
and a fixed topology (i) is given by 

M il \N s ) = D®(N S , iV r )n (i) 

where D^(N S , N T ) is the degeneracy of cluster of size N s , and N T is the 
total number of sites. The total multiplicity of clusters of size N s reads 

M(N S ) = J2 M(i) ( N s) 
(0 

Since the number of spins with a given sign is not fixed in the model (see 
discussion below) one introduces a normalization factor N and 

M(N S ) = ^M(Ns) 
such that 

N T 

N T = N ^ N S M(N S ) 

N S =1 

The interest of this approach lies in the fact that it is easy to work out the 
evolution of quantities like the average size of clusters as a function of the 
temperature 

N T N T 

< N S (T) >= Y, N S M(N S ,N T )/ J2 M(N S ,N T ) 

N s =l N S =1 

the corresponding variance 

a 2 (T) =< N%(T) > - < N S (T) > 2 
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and the fragment size distributions. The expected qualitative features ob- 
tained in a finite system when the temperature of the system increases is 
observed. The results can be interpreted in terms of site percolation con- 
cepts. "Occupied" sites in the percolation framework correspond to cluster 
sites in the present case and "empty" sites to perimeter sites. The relative 
amount of "occupied" sites is governed by the temperature. The number of 
perimeter sites increases with increasing temperature and the relative amount 
of large clusters diminishes. In the percolation picture the number of empty 
sites increases with decreasing occupation probability p, leading to the same 
trend. 

The model which deals with short range nearest neighbour interactions 
has been extended to include long range interactions [ |203|| . Numerical studies 



show that the range of the interaction in finite systems does not change the 
qualitative behaviour of multiplicity distributions in finite systems, except 
for the fact that a long range attractive interaction favours the survival of 
large clusters over a large temperature interval. At high temperature both 



distributions collapse onto each other, as it can be seen in Fig. |15[ The pres- 
ence of "critical" events is seen in the behaviour of the second moment of the 
cluster size distribution. Similar results are obtained with 3d systems. The 
a priori general difficulty to relate lattice models with nuclear fragmentation 
will be discussed later. The use of the Ising model presents a further weak- 
ness which stems from the fact that the identification of clusters in terms of 
spins which are aligned in a given direction breaks the 2.2 (±1) symmetry 
which is present by construction and leads to non conservation of the number 
of "particles" . The remedy to this disease consists of the introduction of a 
term governed by a Lagrange multiplier. We shall come to this point below. 
The present models were thought to build up a bridge between the purely 
geometrical concepts introduced by percolation models and sophisticated 
liquid-gas descriptions. They triggered the construction of more refined mod- 
els aimed to describe and understand the physics of fragmentation and critical 
behaviour in finite nuclei. 

5.2.2 The grand canonical lattice gas model (LGM) 

The preceeding study was aimed to look for the connection between the 
geometrical and probabilistic concepts which define percolation models and 
the thermodynamic properties of spin models which are defined in terms of 
classical Hamiltonians. It did however not explicitly define the link between 
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bond probabilities, energy and temperature. 

This link was established by Jicai Pan and Das Gupta [|205| - |207|| and 
Campi, Krivine and collaborators |208| - |210] . In ref . ||208|| the authors consider 
a classical lattice gas model defined by the Hamiltonian 

N 2 

H — > rij — e > rijUi 

1=1 <«,i> 

where n« = 0, 1 is the occupation number of the lattice site i, N the total 
number of lattice sites and e the energy between nearest neighbours < i,j >. 
The number of particles A = ^\ rii is conserved in the average by means of 
a constraint induced by a Lagrange multiplier procedure. 

The model is treated by means of Metropolis Monte Carlo simulations, in 
the framework of the grand canonical ensemble which can be formally related 
to the canonical description of the Ising model with magnetic field ||2 1 1|| . The 



exploration of the properties of the system in the density-temperature plane 



(p, T) leads to the typical phase diagram shown in Fig. [16| One observes 
the presence of a full line which separates two phases above the line, a gas 
at low p and a liquid at large p. Below the separation line gas and liquid 
coexist in an inhomogeneous phase forming domains (clusters) of different 
sizes. In the thermodynamic limit the transition line corresponds to a first 
order transition. At p = 1/2 in units of the normal density po there appears 
a critical point at T = T c which corresponds to a continuous transition 
characterized by critical exponents which are universal, i.e. do not depend 
on the detail of the interaction between particles. 

The main point of interest in the study of this model lies in the analysis of 
the cluster content of the system for different values of p and T. The identifi- 
cation of clusters raises the question of their definition. One gets confronted 
with two problems [[208|| . The most disturbing one concerns the fact that the 



thermodynamic critical point does not coincide with the percolation critical 
point at which the cluster size distribution shows a power law behaviour if 
clusters are defined as an ensemble of particles occupying connex (nearest 
neighbour) sites ||212j1 . If clusters are defined in a purely geometric way as 
made of connected occupied sites percolation events generate a whole line 
of points in the phase diagram, an infinite cluster is obtained in the infinite 
limit on the right hand side of the dot-dashed line shown in Fig. pl| This is 
of course unphysical, since at high temperature one expects essentially small 
clusters and isolated particles. A remedy to this problem was proposed by 
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Coniglio and Klein [12 13 1 who introduced the energy concept in order to de- 



fine clusters as made of connected ensembles of nearest neighbour particles 
with a bond probability smaller or equal than 

p CK (e/T) = 1 - exp(-e/2T) (30) 

If one uses this definition the thermodynamic and percolation critical point 
coincide and there exists a whole continuous line of points in the (p, T) plane 
called the Kertesz line ||214||, along which the cluster size distribution shows 



a power law behaviour. The important point concerns the fact that the 
definition (|3Up is shown to be very close to the bond probability prescription 
between neighbour particles given by 



p b ( - ) = Prob 



P 2 
2m 



1_ / duu 2 e- u2 (31) 



7T J./e/T 



which stipulates that a pair of particles is bound if (p 2 /2m — e) < where 
the momentum p is drawn from a Maxwell distribution [J20BLEU3]. It is shown 



208|| that if fl3T|) is chosen as a bond criterion then the separation energy S(A) 



of one particle from a bound cluster of A particles satisfies approximately the 
inequality 

S{A) =B{A-1)-B(A) <0 

where B(i) is the binding energy of a system of i bound particles, and the 
"critical" percolation line coincides approximately with the Kersetz line de- 
fined above. 

Following this prescription one obtains cluster size distributions in dif- 
ferent part of the (< p >,T) plane, where < p > is the average density in 
the grand canonical ensemble. These distributions show the same properties 
as those obtained in the framework of percolation models, i.e. a power law 
behaviour corrected for finite size effects along the Kersetz line, with a power 
law exponent r ~ 2.2 — 2.3, an exponential decay at low density and high 
temperature, the characteristic large cluster peak at high density. One may 
notice also that contrary to the case of multifragmentation models J7U], |110[ 



in which clusters are by construction spherical the cluster along the Kersetz 
line have surfaces with fractal properties which allows for a freeze-out den- 
sity larger than the "critical" density. Furthermore one observes that the 
excitation energy per particle above the ground state keeps approximately 
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constant along the critical line | 210|| . Further work relying on lattice models 



in the framework of nuclear fragmentation was developed in the recent past. 
It concerns the study of the two types of phase transitions (thermodynamic 
and percolation) mentioned above which may be of different orders (first and 
second), the importance of the Coulomb interaction, problems related to the 
finite nature of the systems. We examine different approaches and their con- 
tribution to the understanding and explanation of these questions in the next 
subsections. 

5.2.3 Lattice models and their application in the framework of 
nuclear fragmentation 

The generic character of lattice gas models (LGM), their kinship with 
percolation models and classical systems of particles with short range inter- 
actions lead to a large amount of theoretical investigations. 

Das Gupta and Jicai Pan proposed different variants of this type of models 
[|205| , |206| and applied them for direct comparison with experimental data 



||207|| . They worked out the canonical partition function Z on the lattice in an 
analytic form by using the Bragg- Williams approximation for the calculation 
of the potential contribution to the free energy. In this approximation the 
number of sites which are occupied in the neighbourhood of a fixed ocupied 
site is fixed in the average as zn/N, where z is the number of neighbours 
(coordinance) and n/N the average number of occupied sites {n) in a lattice 
with iV sites. The approximate analytical expression of Z allows to determine 
different equations of state, in particular the relation between the pressure 



and the volume which comes out as a van der Waals-type equation |205 



In a further development [|206|| the authors introduced a grand canonical 



description, the mean field Bragg- Williams approximation was replaced by 
the Bethe-Peierls approximation which corresponds to the introduction of 
blocks of lattice sites inside which occupied nearest neighbour sites interact. 
An interaction e between blocks is taken into account and fixed along with the 
chemical potential of the particles by means of a self-consistent procedure. 
This further improvement brings in fact rather small changes when compared 
with the Bragg- Williams approximation and is rather close to the outcome 
of a mean field approach using a realistic Skyrme interaction. In all cases the 
descriptions show signs for the existence of a first order phase transition for 
all values of p, a second order transition at p c = 0.5po where po is the normal 
nuclear matter density. 
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Emphasis was put on the study of the cluster content of finite lattice gas 
systems. It was shown by means of explicit calculations that the second mo- 
ment of the fragment size distribution shows a maximum at different values 
of the temperature for different values of the density |205|| . In an interpre- 



tation in terms of percolation concepts this is the sign for the existence of a 
Kertesz line. If one introduces 

P 2 r /2fi-e < 
or > 

(where p r is the relative momentum between neighbouring particles, \x the 
reduced mass and e the strength of the interaction) as a criterion to decide 
whether the considered particles are bound or not the bond probability is 
the same as the one given by (|3l"D in section 5.2.2 ||206|| . Effective exponents 



r extracted from the experiment in the spirit of refs. [^T|,^] (see section 2.1) 
were compared with lattice calculations in an attempt to obtain informations 
about characteristic freeze-out densities of systems formed at different beam 
energies. The freeze-out density determined in this way came out to be 
P ^ 0.4p . 

One very important point in nuclear fragmentation concerns the role and 
importance of the Coulomb interaction. This point has been emphasized 



in the framework of multifragmentation models [7JJ. A calculation in the 
framework of an LGM for a large system obtained through the collision of 
Au on Au shows that the Coulomb interaction may have a sizable effect 



and may lead to large discrepancies between LGM and the experiment ||207 
One may introduce two types of binding energies between nearest neighbours 
which correspond to the interaction between identical and different types of 



particles (protons and neutrons) ||215|| . This has been done by means of 



Monte Carlo simulations and the results were compared to the outcome of 
molecular dynamics (MD) calculations which coincide with the lattice gas 
results in the absence of the Coulomb interaction. It is no longer the case 
when this interaction is taken into account. The LGM results deviate from 
the MD results in large systems with a sizable amount of protons. We shall 
come back to the Coulomb problem below. 

Further comparisons with different models have been performed | |216| | such 



as percolation and multifragmentation models. The analysis of different ob- 
servables shows that these models share common features, in particular the 
statistical multifragmentation and lattice gas approaches, except for quanti- 
ties like the specific heat Cy at low temperature. The multifragmentation 
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approach possesses quantum mechanical properties and Cy increases linearly 
with T while in the LGM Cy stays at the constant value 3/2. It should be 
mentioned that the algorithm used in the simulations |[205| - |207| does not 



sample the events in a rigorous fashion [[217|| . This should however not qual- 



itatively affect the results and change the conclusions. 

5.2.4 Extensions of the LGM in the framework of the canonical 
ensemble 

Lattice gas models were further developed and analysed in the frame- 
work of the canonical ensemble by Gulminelli and Chomaz ||218| - [221[| , Borg, 



Mishustin and Bondorf ||222|| and Carmona, Richert and Tarancon [|197|| . 



Gulminelli and Chomaz | |218| , |220| constructed the phase diagram of a 
finite system by working out the compressibility coefficient. This quantity 
shows an anomalous behaviour, it becomes negative at some temperature T 
over some density interval, the sign for the existence of a first order phase 



transition in the infinite system. Fig. [17| shows the behaviour of the chemical 



potential represented as a function of pressure exemplifies this phenomenon. 
It has also been observed that the fragment size distribution possesses a 
separation line (Kertesz line) which crosses the phase separation line and 
enters the coexistence zone at some low value of the density. This is a priori 
somewhat surprising since it means that in the thermodynamic limit one 
would observe the existence of a second order transition in the coexistence 
zone. However, it is found out that for large systems the separation line stops 
effectively at the thermodynamic critical point located at p/po = 0.5, hence 
the authors interpret the effect as a typical finite size effect. They comfort 
this point of view by showing that in this region of the phase diagram the 
fragment size distribution does not go through a power law distribution when 
one crosses the separation line in the coexistence zone. If this is so, it means 
that the system never shows the signs for a critical behaviour for subcritical 
densities. 

Lattice gas models were also recently extended to the case where two 
types of particles, protons (p) and neutrons (n) coexist. The refinement is 
implemented by the choice of different nearest neighbour interaction strngths, 
£np 7^ £nn = £p P - These constant quantities were chosen such as to reproduce 
the nuclear binding energies of the corresponding quantum systems. One may 
ask for the consequences of the introduction of the isospin degree of freedom 



||219|| . In fact, the extension of the model has only a small effect on the phase 
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diagram, so that in practice there exists only one relevant order parameter, 
Pl~ PGi where pi, and pc are the liquid and gas density respectively Critical 
exponents do not depend on the isospin degree of freedom either, and the 
critical temperature shows a quadratic dependence on asymmetry. The sign 
for isospin effects may perhaps be found in the fact that in the coexistence 
region the vapour mode of small fragments is more isospin asymmetric than 
the liquid mode of larger fragments. This can in fact be observed in the 



experiment 268.269 



A similar model with protons and neutrons has been studied in ref. ||222 
The thermodynamic properties of the system are determined for open bound- 
ary conditions and the caloric curve is constructed at constant density and 
constant pressure. For fixed pressure one observes clearly the appearance of 
a plateau which is the sign for a first order transition. For fixed density, the 
caloric curve is monotonously increasing with energy. Fragment size distri- 
butions are analysed along isobars and isochores. They show the usual trend, 
with increasing temperature the large fragment yield shrinks and the small 
fragment and light particle contribution increases continuously. 

Recently the effect of the shape of the surface of the fragmenting system 
on the fragment size distribution has also been investigated [ 221|| . Cubic 



shapes with and without periodic boundary conditions as well as spherical 
shapes with sharp and diffuse surfaces have been considered. These shapes 
have no real influence on the thermodynamic properties and fragment size 
distributions of small systems with typically 100 particles at the critical den- 



sity. Similar tests have been performed on other models, see below ||224 
It raises the question wether the shape dependence of fragmentation events 
may be detected in experimental data. 

5.3 Cellular models of nuclear fragmentation 

As already mentioned the success of percolation concepts for the repro- 
duction of experimental fragment size distributions raised the question of 
the correlations between theses quantities and the actual dynamics of highly 
excited nuclei |169| , |203(| . Peripheral collisions |225|| indicate that the spec- 



tator remnant behaves essentially as an excited system of particles which is 
thermalized. It leads to the idea that fragmenting nuclei can essentially be 
described as completely disordered systems of interacting particles, at least 
under the experimental conditions described above. Furthermore, it is a non 
trivial but established fact that disorder may be connected with complexity 
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and universality ||226||. This property reflects the fact that disordered systems 



are not sensitive to the details of their dynamical characteristics but show 
features which are common to many of them such as the scaling properties 
discussed above. 

In the present section we introduce and discuss so called cellular models 
which are aimed to relate the phase space properties of a classical system of 
particles to bond probabilities and hence to the generation of bound clusters. 
They can be considered as extensions of lattice models. 

5.3.1 Disordered systems and cluster identification 

A fragmenting nucleus is considered as an assembly of interacting particles 
contained in a finite three dimensional volume V divided into A cubic cells, 
V = Ad 3 where d is the linear size of each cell. Each cell contains one 
particle | |224| , ^27|] . Each particle inside its cell is characterized by its position 
Tiii = 1, . . . , A) and linear momentum pi, both phase space coordinates 
being taken at random and obeying \fi — Tj\\pi — Pj\ > h which has to be 
verified for particles located in neighbouring cells. The interaction is chosen 
as a short range two-body potential adapted to a classical description of 
the nucleus | 228| , 229 |. It acts between nearest neighbours only. The total 



energy of two such particles can be split into a centre of mass and a relative 
contribution 

E = P 2 /2M + p 2 /2fi + V(r) = P 2 /2M + e 

with P = Pi+ Pj, p = Pi — Pj, r = \fi — fj | , M = 2m and /i = m/2, where m 
is the mass of the nucleon. Each coordinate r*j is defined as Xi = x i0 ± r/d/2 
(same for %ji and Zi), where x i0 locates the centre of the cell and r\ is drawn 
uniformly in the interval [0, 1]. The momenta {p{\ are taken uniformly in a 
sphere of radius p ma x- Two neighbouring particles are said to be bound to 
each other if e < and unbound if e > 0. By means of numerical simulations 
the criterion allows to construct a bond probability p(E) = n(E)/'N(E) where 
yf(E) is the total number of draws in an energy interval [E, E + 5E] and n(E) 
the total number of events out of N(E) for which the two-particle system is 
bound at energy E. The violation of the inequality relating relative distances 
and momenta defined above leads to the rejection of the corresponding event. 
Three different procedures have been introduced in order to identify 
bound clusters. In the first called A2, the energy E of each pair is drawn from 
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the distribution 'N(E) and the corresponding p(E) is compared to a random 
number rj belonging to [0, 1]. If p(E) > r\ the pair is considered as bound 
and unbound in the reverse case. This procedure is close to a percolation 
algorithm but p[E) is not chosen in a uniform distribution. It neglects also 
the correlations which exist between pairs sharing common particles. 

The second procedure (Al) fixes the energy of a pair E = E{fi,Pi, Tj,Pj) 
of neighbouring particles % and j by means of their phase space coordi- 
nates which are fixed as described above. This fixes also p(E) and the 
bond criterion is the same as for A2. Finally, correlations are fully taken 
into account by means of algorithm A0 for which the generated phase space 
{fi,Pi,i = 1, . . . , A} is used to calculate explicitly e for each pair. It cor- 
responds to a deterministic determination of the clusters for each generated 
event. 

It is important to notice that none of these descriptions implies that the 
considered systems are necessarily in any kind of equilibrium. 

5.3.2 Physical observables. Confrontation with bond percolation 

The fragment content can be characterized by the mass distribution, its 
moments and associated observables like 72 which has already been defined 
above. In order to compare directly with the experiment these observables 
are determined as a function of the multiplicity to which can be measured 
and is univoquely related to the bond probability p. For each event k 

i 

or S 2 (to) = m 2 (m)/m{ (to) 

where P^ (to) is the number of clusters of size i and the prime index indicates 
that the heaviest fragment is omitted. Similarly 

7 2 = m (m) ■ m 2 [m)/m\ (m) 

= 1 + cr ' (m) • m (rnj/mi (m) 

where a^' = I raj, (to) • to 2 (to) — m\ (to) i/toq (to) 



Average quantities are defined as 

< m 2 ( m ) > = N' 1 J2 



Nev 

m 2 (m) 



k=l 

< 72 (m)> = N-JJ2^2 k) 



m) 
fc=i 



Numerical investigations have first been performed with the algorithm A2. 



The comparison with bond percolation is shown in Fig. [Tj|. One observes an 
amazing agreement between the two descriptions although bond probabilities 
are taken from different distributions and this agreement is very robust with 
respect to the strength of the two-body potential. In fact, if p{rj) is the 
probability distribution of r\ = p(E ) for a given E the equivalent probability 
in bond percolation is 



rv 
fj — 1 — / dr]'p(r]') 
Jo 



Through this relation m and f\ are related to each other in the same way as 
m and p in bond percolation. In particular, one observes a critical value fj c 
for which an infinite cluster appears in the infinite system. 

The agreement survives also if protons and neutrons are differentiated 
by means of different bond probabilities p pp (E) = p nn (E) ^ p pn (E) |224| |. 
Qualitatively similar but quantitatively different results are obtained with 
algorithm Al, in particular the location of the maxima of < ui2 > and 
< 72 > are shifted to higher values of m. The deviations are indicated by 
the correlations which are built in Al and absent from A2. One expects 
that the effect is due to the hybrid nature of the algorithm which mixes 
probabilistic and deterministic properties. 

Algorithm AO provides a priori the most realistic description since it takes 
all correlations consistently into account. It allows also for an easy treatment 
of the long range Coulomb interaction. For a fixed value of the length d of the 
cubic cells the moments < m 2 > and 72 =< m 2 > • < m > / < mi > 2 can 
differ sensibly from the percolation results. If, however, one averages over 
events corresponding to values of d which are uniformly distributed over a 
given interval of densities p, say from 5 x 10~ 3 to 0.17 fm~ 3 , one obtains mass 
distributions, < m 2 > and 72 which are remarkably close to those obtained 
with A2 and robust with respect to modifications of the strength of the short 
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range potential as it is the case with A2. The agreement between the results 
obtained with AO and A2 is not easy to understand. It seems to indicate 
that averaging over densities acts as a mean to wash out correlation effects. 

5.3.3 Confrontation with experimental data 

Multifragmentation data available from experiments made by the 
ALADIN collaboration [178, 230| , 231 1 have been confronted with the out- 



come of the A2 algorithm in the same way as it was done for the percolation 



approach in section 4.7 |224| , |227| . The charge Z to mass A relation used is 
the same as in this section. For fragments of mass A = 3 the isobaric ratio 
3 H/ 3 He is fixed to 2.2. The mass of the spectator part of the fragmented 
system depends on the impact parameter, the number of particles involved 
is divided into 10 different sets corresponding to 10 intervals for the impact 
parameter. Each system with Ai (i = 1, . . . , 10) particles is put in a cube of 
7 3 cells in as compact and symmetric a geometry as possible. Each calculated 
observable O is averaged over the whole impact parameter range 

10 / 10 



< O >= J2^Ab t O(bA / J>A&; 



Results of this analysis are shown in Fig. [19] for different observables, the 
average value of the charge of the largest fragment Z max , the IMF (3 < 
Z < 30) multiplicity, the He isotope multiplicity which are represented as a 
function of Zf) 0un d and the average number of fragments with charge 6 < Z < 
30. As it can be seen, the agreement between experiment and simulations is 
excellent. This is not really surprising since bond percolation and A2 lead to 
very close results and, as it has been seen in section 4, percolation reproduces 
experimental results. 

5.3.4 Geometric and collective effects 

The use of simple cubic shapes for the description of the fragmenting 
system may appear as too restrictive. This point has been investigated by 
introducing other shapes, spheres, parallelepipeds, streaked, crenellated and 
crumpled surfaces which were generated by means of a random growing pro- 
cess leading to fractal surface structures which may introduce empty cells in 



the bulk of the system [|224j| . The moments of the mass distributions show 

80 



deviations from the cubic case, but these deviations are always surprisingly 
small, they never lead to qualitative changes in any observable related to 
fragment size distributions. 

Cellular models describe disordered systems which show a power law be- 
haviour with deviations due to finite size effects at a value of the fragment 
multiplicity corresponding to criticality in the infinite system. Experimen- 
tally this situation prevails for peripheral collisions ||178| , p30| , |231|| . In the case 



of more central collisions fragment size distributions show rather an exponen- 



tial fall off as a function of the mass number ||232|| . These collisions are also 
characterized by a sizable collective radial flow |155| , [233f . It is out of scope 
to reproduce this dynamic effect with the present simple-minded approach. 
It is however possible to mimic it. This was done with algorithm AO, mixing 
events with linear cell size 1.8 < d < 6 fm for A = 216 particles. The linear 
momenta strengths of the particles are chosen randomly but their directions 
are fixed along a line which joins the centre of the system to their actual 
position, pointing out of the system. The typical behaviour of the mass dis- 



tribution in this case is shown in Fig. 20, along with the case where collective 



effects are absent. The shape which is power law like in the case where the 
directions of the momenta are random gets very close to an exponential for 
sizes between 20 and 90 units of mass. Deviations are observed for light and 
heavy fragments. Such a behaviour has been confirmed in ref. ||155| were a 



more detailed analysis on a more sophisticated approach has been performed. 
This result comforts the idea that peripheral collisions induce the generation 
of full disorder and are very weakly affected by collective effects whereas 
central collisions show sizable collective effects. 

5.3.5 Thermodynamic properties of systems described by cellular 
models : an analytical model 

The encouraging results concerning fragment multiplicities obtained in 
this type of simple models raises the question whether the present description 
supports the existence of a thermodynamic phase transition as it is the case 
for the lattice gas models described above. The problem which comes up is 
the fact that the generated events do not need to possess the property of being 
in thermodynamic equilibrium since the algorithm which is used does not 
guarantee it. Indeed, a necessary condition for thermodynamic equilibrium 
in the canonical ensemble is the determination of an absolute minimum in the 
energy for a fixed temperature through the application of detailed balance. 



This has not been implemented. However there are reasons to believe that the 
generated systems are close to if not at equilibrium. This led to the analysis 
of their thermodynamic properties in the framework of the microcanonical 
ensemble ||234|| . The entropy S(E) has been determined as a function of the 
energy E and the caloric curve constructed by defining the temperature as 
T _1 = (dS/dE)y where V is the volume of the system. The outcome of the 
study is the fact that T increases monotonously with E and there is no sign 
for the presence of a transition. An analytical investigation which describes 
the system in a way which is close to the one which is obtained through 
numerical simulations confirms this result [ |234| . 

The absence of a transition is in fact qualitatively understandable if one 
goes back to the lattice gas models of section 5.2. Indeed, the discussed 
cellular models correspond to the extreme case where the system is homoge- 
neously filled with particles, hence p/po = 1 in Fig. [TJ] where the system is 
in a pure liquid phase. 

In order to study the behaviour of a system for every density one may 
introduce a simple 3 dimensional cubic lattice whose volume V is divided 
into N 3 cells of linear size d, V — (Nd) 3 . This geometry can be straightfor- 
wardly extended to a parallelepiped V = NiN 2 N 3 d 3 |[235|] . The restriction to 
attractive nearest-neighbour interactions leads to a Hamiltonian 

i=l <n.n>i 

where A is the total number of particles, < n.n >k stands for the "fcth 
nearest-neighbour" sites(i, j), Vq is the strength of the interaction and Sk = 
(resp. 1) corresponds to a cell which is empty (resp. occupied by a particle). 
The model which describes an inhomogeneous system can be trivially mapped 
into a spin-model by means of the change of variables 

o"fc = 2sfc — 1 
Particle number conservation requires 

N 3 

1/2^(1 + ^) = A 

fc=i 
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In terms of the new variables the Hamiltonian reads 

N 3 



2m 4 _ 

fc=l <n.n>i fc=l 



*° = E&-t[£ «^+«E*+^ 



which corresponds to an Ising model with fixed magnetic field h = 3Vo/2. In 
order to work out its thermodynamic properties one goes over to the contin- 
uum limit by replacing o k = ±1 by the continuous variables u k in [— oo, +oo]. 
This leads to a description in the framework of the spherical model which 
is exactly solvable in any number of dimensions for different types of inter- 
actions and even in the presence of disorder [ |236| - |240f . Constraining the 
continuum variables {u k } by means of 

N 3 

J2 u l = N * ( 33 ) 

fc=i 

and introducing the Lagrange multipliers A and /i leads to the Hamiltonian 

N 3 N 3 

fc=i fe=i 

and a canonical partition function 

/+oo _ _ /"+oo 

Y[dp% du k exp[-(3H({u k })] 

■°° k,i J -°° 

where p\ corresponds to the zth component of the moments p k of the particle 
k. The parameters A and /x are fixed by 

- P^dhiZ/dX = N 3 (34) 

-P~ x d\D.Z/dn = A (35) 

The physical justification for the introduction of the continuum limit variable 
{u^ is the fact that the interaction in Hq can now be interpreted as taking a 
continuum of values, not only ±Vo/4, which is sensible in the case of cellular 
models, since the location of the particles in cells is not restricted to the centre 
of the cells. Because of the constraints imposed by the volume constraint (|33|) 
the important contributions to V <n , n>1 remain finite. The procedure does 
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not take care of the fact that positions of a set of interacting particles are 
correlated which induces correlations on the two-body contributions. These 
correlations are not taken into account here, like in the use of algorithm A2 
in numerical simulations. 

The partition function can be worked out explicitly. It can be cast in the 
form 

Z{/3, A, N 3 ) cc exp[-pN 3 f(P, A, N 3 )} 

with the free energy per cell 

N 3 -l 

f(J3, A, N 3 ) = [(/i - 3V )/2] 2 /[A(3V /2 - A] + (1/2(3N 3 ) ■ ]T ln[/3(A - V a t /2)] 



i=0 



where a* = «ni 2 i 3 = }] cos & fc 

fe=i 
& fc = 2m k /N {i k = 0, 1, . . . , iV - 1) 

The constraints d34T) and (1351) lead to the relation 



N 3 -l 

J2 [A - W2]" 1 = 2iV 3 /5(l - M 2 ) (36) 

i=0 

which fixes A. Then 

H = 4(3V /2 - \)M + 3V (37) 

where 

M = -[l + 2((3N 3 )- 1 d\nZ/d/j] 
= -0*-3F o )/[4(A-3Vo/2)] 

= -l + 2A/iV 3 

= 2p/p -l 

which shows that M G [—1, +1]. 

The dispersion relation (^TJ) has iV 3 poles. There exists one solution 
which corresponds to A > Voa ma x/2 with a max = 3 and guarantees that the 
corresponding value of A leads to an expression of Z which makes physical 
sense. 
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It can be shown [239 that for M = (p/po) — 0.5) the free energy per 



unit volume gets singular for a given value of the temperature T = T c (0) 



V /2W 3 (0) where W 3 (0) = 0.2527 E3a,E33|. The same behaviour of the free 



energy is observed for all other values of M leading to critical temperatures 
T C (M). The specific heat is given by 

C v = -(3 2 de/df3 

where e = E/A is the energy per particle. Cy shows a maximum for every 
value of M at T C (M) which signals a continuous phase transition between two 
phases. The corresponding phase diagram is shown in Fig. [Hj. The analytic 
expression of Z allows for the determination of all thermodynamic quantities 
and the construction of the caloric curve T = f(E) which is a monotonously 
increasing function. This is due to the fact that Cy remains finite at the 
transition point for any value of p/p ||235|| . Only the derivative of Cy shows 
a discontinuity at the transition points. 

The model can be extended to the case where particles interact through 
short or/and long range interactions 



H = K - [V 01 y^ SiSj + V 02 ^ s 'i s i H h V 0m 5Z SiS > 

<n.n>i <n.n>2 



<n.n>i <n.n>2 <n.n> r 

2 



where K is the kinetic energy and < n.n >& has been defined below (p2|). A 
sensible parametrization in the case of nucleus would correspond to 

V&2 = V^ 



^Om — V { 

(k) 

where V s > stands for the short range nuclear potential and V l < for 
the long range repulsive Coulomb interaction parametrized in the following 
form 

Vt k) = (Z p /A)e 2 /kd 

where Z p is the total number of protons. 
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The partition function can be worked out in the framework of the spher- 
ical model in the same way as in the former case |[235|| . Fig. [21] shows typical 



phase diagrams for different cases involving long range interactions. As one 
can see, they lead to qualitatively similar results. 



In the usual spherical model in which the constraint fl35|) is not imple- 
mented p plays the role of a fixed magnetic field. Then, for /i / 0, the 
magnetization can jump from — M to +M which corresponds to a first order 
transition. For p = the transition is continuous. This is not the case here 
since the magnetization M (density p) is fixed, and not p which is determined 
through fl3~5|). It explains why the transition is continuous in the present case 
for any value of M or (p/p ). 

5.3.6 Thermodynamic properties of systems described by cellular 
models : numerical simulations 

The former approach is a simplified even though quite realistic thermo- 
dynamic description of a disordered system. However, it does not allow to 
work out the fragment content of the system, which are accessible by means 
of numerical simulations ||241|| . This has been performed in the framework 



of the canonical ensemble for particles which experience short range nearest- 
neighbour interactions. 

The Hamiltonian is written as 



A 2 



2m 

=1 <i-j>i 



where Vu is the potential which acts between particles located in neighbour 
cells ||228|| at a distance r^ from each other, A is the number of particles, p 



A/L 3 the density of particles in a volume V = L 3 and the s^'s are occupation 
numbers defined above. Periodic boundary conditions are imposed in the 
numerical simulations. 

A Metropolis Monte Carlo procedure is used in order to generate events 
corresponding to different densities p and temperatures T. For each (p, T) 
the total number of events is 6 • 10 4 and for each event about 5000 Metropolis 
iterations have been realized. Each iteration consists of either the move of 
a particle from its original cell to an empty cell, or its displacement inside a 
cell, or the determination of a new momentum. Each operation is realized 
with a probability 1/3. 
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The latent heat is defined as 

C v = j(<U 2 >-<U> 2 ) 

where the average is an ensemble average and U the potential energy of the 
system. 

The calculations have also been extended to the microcanonical ensemble 



with fixed A and energy E, by means of an algorithm described in ref. [ [242 
which is well adapted to the present type of models. The probability for the 
existence of an event which has a configuration s characterizing the positions 
and momenta of the particle is given by 



av -i 



W E (s) oc(E- U{s)) 2 Q(E - 17(a)) 

where d is the space dimension, V the volume and the Heaviside function. 
Events are selected by means of a Metropolis procedure which accepts a 
configuration s' starting from a configuration s with a probability 

W x -+ X > = mm 1, 



W E (s) 
The temperature is obtained from 

T = 2 < K > /Vd 

where < K > is the average kinetic energy of the particles. 

Typical examples of the behaviour of Cy and the caloric curve are shown 
in Fig. ^. The specific heat presents a more or less pronounced enhancement 
for some value of T. This enhancement is taken as the finite size restriction 
of a point which determines the phase separation line between a liquid and a 
gas in the thermodynamic limit. The numerical simulations show a smooth 
monotonous increase of the energy with temperature for fixed volume which 
is very similar to the results obtained in the framework of liquid gas models. 
For these small values of A and V it is not possible to identify the order of 
the transition. It is also remarkable to see that canonical and microcanonical 
calculations are very close to each other, even though the size of the system 
and the volume which contains it are small. We shall come back to this point 
and develop it in section 6. 

The identification of bound clusters leads to the separation of the frag- 
ment content into two areas of the (p, T) phase diagram separated by the 
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equivalent of a Kertesz line, above which the system contains essentially in- 
dividual particles and small fragments, whereas large clusters are present 
below it. The separation line is expected to correspond to a power law frag- 
ment mass distribution in the limit of an infinite system. The presence of 
the power law behaviour is detected by analysing the calculated mass distri- 
butions. 

The results are shown in the phase diagram of Fig. |23|. Several remarks 
are in order. The thermodynamic separation line shows a dip at p = 0.5. This 
is a finite size effect. A study of this point on a spin model which presents the 
same property confirms this point, see below in section 6. The separation 
line itself is not rigorously determined, there exists other means to fix it 
and furthermore it does not correspond to the limit of the infinite system. 
Contrary to the result of ref. |[218|| , the Kertesz line stays above the separation 
line for all values of p but does not end at p — 0.5. Furthermore, at p = 0.5 
the thermodynamic critical point does not coincide with the corresponding 
point on the Kertesz line. The reasons for this are not easy to identify clearly. 
It makes sense to believe that all these points are related to finite size effects 
and possibly also the numerical precision with which the relevant quantities 
can be obtained. The determination of critical exponents by means of finite 
size scaling arguments has not been performed. The analysis and outcome 
would be similar to the ones obtained in the framework of another model 
which will be discussed below. 

The remarkable point which may be retained from this study is the fact 
that the a priori more sophisticated cellular models discussed here are very 
close to the more schematic lattice gas models. If there are quantitative 
differences it is, at present, out of scope to judge which is quantitatively 
closer to the experiment. 



6 Finite nuclei and phase transitions : 
present experimental and theoretical sta- 
tus 

The experimental quest for a thermodynamic phase transition in nuclear 
matter received a new impetus in 1995 with the work of Pochodzalla et al. 
who tried to construct the so called caloric curve 112431 . If nuclear matter 



would behave like a Fermi gas one would expect an increase of the temper- 
ature T with excitation energy E* such that T oc E* 1 ^ 2 . Deviations could 
indicate the existence of a phase transition which, following classical argu- 
ments, is expected to be of the liquid-gas type. 

However, nuclei are finite and the asymptotic limit of infinite matter 
is experimentally out of reach. The situation which is common to other 
fields like physics of aggregates and mesoscopic systems in condensed matter 
physics needs the development of new tools which may allow to detect the 
existence of the transition through the study of the behaviour of physical 
observables which characterize finite systems. 

In the present section we try to give an overview of the recent experimen- 
tal efforts which have been made and the theoretical investigations which 
have been initiated, specifically with the help of lattice-type models in order 
to detect a possible phase transition and unravel its nature in excited nu- 
clear matter under the assumption that experimental measures are made on 
thermodynamically equilibrated systems. 

6.1 Caloric curves and phase transitions : experimen- 
tal status 

The caloric curve which relates the internal energy of an excited system 
at thermodynamic equilibrium to its temperature is a priori the simplest 
experimental tool to look for the existence of a phase transition. The first 
determination attempted by the ALADIN collaboration |j243|j collected the 



outcome of 197 Au on 197 Au collisions at 600 MeV-A with data obtained by 
means of less energetic collisions. The excitation energy was determined by 



following the procedure prescribed in ref. |182|| and the corresponding tem- 
perature fixed by means of arguments relating this quantity to the so called 
double ratio procedure |244j| , in the present case the ratios of 3 He/ A He and 



& Li/ 7 Li isotopes. The corresponding curve showed the features of a rather 

89 



strong first order transition, with a characteristic close to constant temper- 
ature T over a large energy interval lying between 3 and 10 MeV excitation 
energy per nucleon, which may be interpreted as a sign for the generation 
of latent heat and followed by a strong increase of T with excitation energy 
above 10 MeV. A critical discussion followed this observation, in which the 
hypotheses and simplifications underlying the definition of the temperature, 
the freeze-out density and the increase of T at high excitation energy were 
examined |245|| . 



A similar analysis by means of 1 GeV-A collisions of C on Au was initiated 
soon after this first attempt ||246f| . The reaction was interpreted as a two- 



step process, a prompt stage with the emission of light particles and a second 
step corresponding to the decay of an equilibrated system undergoing frag- 
mentation. The outcome of the analysis lead to a smooth caloric curve from 
which any clear sign corresponding to a first order transition was seemingly 
absent, showing at most the possible existence of a continuous transition. 
In this experiment the temperature was determined in terms of two different 
isotope ratios, ( 3 ife/ 4 i7e)/( 6 Li/ 7 Li) and ( 3 He/ 4 He)/(d/t). The results were 
essentially confirmed later in a further experiment ||247|| . 



The INDRA collaboration considered Ar on Ni collisions at 52 and 



92 MeV ■ A in order to investigate this point ||248|| . Several double isotopic 
yield ratios were used in order to define the temperature. They led to differ- 
ent slopes in the temperature increase as a function of the excitation energy. 
The issue of the experiment and its interpretation have been performed by 
means of statistical models [ |138| , |139| , |142| . A criticism has however been 



raised |249| which claims that the observation of the transition may have 
been missed because of the choice of too large steps in energy binnings. 

At this stage the ALADIN collaboration considered the problem of the 
determination of the temperature by means of an analysis of Au on Au colli- 
sions at 1000 MeV -A p25|| which was determined through different isotopic 



yield ratios and corrected for sequential evaporation from primary fragments. 
In this experiment, temperatures were also defined by means of excited state 
populations [|250|,|251|. Both types of temperatures are shown in Fig. 123. 



Those obtained from population yields are remarkably constant as a func- 
tion of the excitation energy per particle, in contrast to the temperatures 
obtained by means of isotopic yield ratios. The interpretation of this appar- 
ent inconsistency leads to the conclusion that these temperatures may cor- 
respond to the final stage of fragment emission. The corrected new caloric 
curve emerging from the use of isotopic yields does no longer show a clear 
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plateau-like behaviour. At the highest value of the excitation energy there 
remains however a sizable rise in the temperature which seems to indicate 
the entrance into a new phase. 

A further attempt has been performed by bombarding Ag and Au with 
4.8 GeV 3 He ions |252fl . The temperature was fixed by means of the dou- 
ble ratio yield ( 2 H/ 3 H)/( 3 He/ A He). A slope discontinuity is observed for 
E*/A ~ 2 — 3 MeV which is followed, with increasing excitation energy, by 
a monotonic rise up to the measured energy E*/A ~ lOMeV. The curve is 
in qualitative agreement with the theoretical predictions of the Expanding- 
Evaporating Source model (EES) 



^, [253| and the Statistical Multifragmen- 

tation Model (SMM) jTT0| , g5i . 

More recently, a further result has been obtained on excited systems with 
A ~ 110 particles produced by means of different projectiles and targets with 
a bombarding energy of 47 MeV ■ A [ |255| |. The temperature was obtained 
from the double ratios (d/t) /{^He/ A He). An analysis based on coalescence 
models [14C, 25^-259] was used in order to follow the evolution of temperature 
and densities. Calorimetry measurements fix the excitation energies. The 



experimental points combined with former results |260|| show a steep increase 
of T up to E*/A ~ 3.5 MeV. The temperature stabilizes for 3.5 < E*/A < 
7 MeV around T ~ 7 MeV. 

The analysis and comparison of results obtained by different groups show 
hardly more than qualitative similarities and are rather non conclusive as far 
as the existence and possible order of a phase transition are concerned. They 
are plagued by different problems which are not only related to finite size 
effects or the equilibration problem, but also by the fact that in the exper- 
iment the system undergoes expansion, sequential decay of primary excited 
fragments, apparently temperature determinations can lead to inconsistent 
results. These problems have been investigated by several groups [ 261 267 . 
Further new methods for the determination of the temperature have been 
proposed |[268 |. As the matter stands, there are signs which point towards 
the existence of a transition. The caloric curve may however not be the most 
adequate tool to confirm it. We shall discuss this point in the next section. 
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6.2 Signs for phase transitions in finite systems : tests 
and applications to the nuclear case by means of 
lattice models 

Nuclei are finite systems, hence experimental observations refer to small 
systems, far from the thermodynamic limit. Contrary to other aggregates or 
quantum dots whose size can be experimentally varied and controlled, there is 
no such flexibility in fixing the number of constituents over large scales here. 
It shows that even if the experimental determination of physical observables 
could be fully tested, there would remain the problem of the small size limit. 
This point raises the question of the detection of the existence and nature of 
a phase transition through the measurements on finite systems. The signs for 
the existence of transitions clearly exist, they are seen in finite systems since 
infinite systems do not exist, and one may speak about phase transitions 



in finite systems p70|| . The determination of their order necessitates good 
tests and a great deal of caution. This point has been raised in several fields 
of physics in the recent past. We review some of the methods and tests 
which have been proposed in the next section before applying them to lattice 
models and the analysis of experimental data. 

6.2.1 Tests for the existence of phase transitions from finite sys- 
tem analysis 

There exist several characteristic properties which may signal that a finite 
system can show a critical behaviour related to a phase transition. 

Close to a critical point observables follow scaling laws which relate their 
behaviour at different scales, i.e. for different sizes of the system. They allow 
for the determination of critical exponents which characterize the singular 
behaviour at the critical point in the limit of an infinite system. This property 



has already been presented and discussed in subsection 4.3 above |197|| . We 
shall come back to it below. 

A second characteristic of finite systems which are close to a phase tran- 
sition is the enhancement and universal behaviour of fluctuations in ob- 
servables and order parameters. This type of properties is commonly stud- 
ied |271|| and has been used | |272| | as will be shown below. Recently the prop- 



erties of order parameters of finite systems close to a second order transition 
have been studied in the framework of a finite size scaling approach ||273||. It 



comes out that order parameters can be characterized by probability densi- 
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ties which possess specific properties depending on their location with respect 
to the critical point. In particular, at the critical point, they show universal 
scaling and deviations from a gaussian behaviour in the tails. This behaviour 
is valid for both equilibrium as well as off-equilibrium processes in systems 
which undergo second order transitions. The method could work as an alter- 
native to the determination of critical exponents. Applications will certainly 
be developed in the near future. 

Informations about the thermodynamic properties of small extensive and 
non-extensive systems (for which the sum of the energies and entropies is not 
the sum of energies and entropies of their parts, f.i. if long range interactions 
are present) can be very efficiently analysed in the framework of the micro- 
canonical ensemble. This aspect has been extensively discussed, see f.i. p7(J| | 
and refs. therein. The important point of practical importance concerns 
the topology of the entropy S(E, A) of a system with A particles and total 
energy E ||270|| . In the presence of a first order transition, the entropy shows 



convexity (negative curvature) in certain ranges of energy. This behaviour 
is specific to the microcanonical ensemble which is able to give a faithful 
description of inhomogeneous systems. In the caloric curve the temperature 
is a multivalued function of the energy, showing a bend ("S"' curve) and, 
as a consequence, the specific heat can get negative over some intervals of 
energy. The physical explanation for this phenomenon is the existence of a 
surface tension which is generated by the coexistence of two phases. Indeed, 
the creation of intra-phase surfaces generates an entropy loss ASsurf which 



can be obtained from the relation |274j| 



A 9 - 7 S 

where 7 is the surface tension, T c the transition temperature, £ the surface 
area of all species with mass larger than 1, and A the number of particles. 
This remarkable property shown by finite systems will be used below as a 
main tool in order to detect and characterize the behaviour of fragmenting 
nuclear systems. 

The present collection of tests is certainly not exhaustive, there may ex- 
ist other means to detect signs for phase transitions. One of them which 
grounds on the Yang and Lee theory has already been mentioned in sec- 
tion 5.1.3 ||199|| . Recently Ma ||275| , |276f proposed to use the enhancement of 



entropy production at the pseudo-transition point and claimed that Zipf's 
law IP 77(1 is verified at the transition point. This has been tested numerically 
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by means of calculations involving an LGM model and Molecular Dynamics 



simulations. Ref. ||278|| discusses the correlation between an "S"' curve in 
the caloric curve in the microcanonical framework, the thermodynamic tem- 
perature T = (dE/dS)v and the temperature determined from the kinetic 
energies of the constituents of the system. It is shown that the caloric curve 
obtained by means of the thermodynamic temperature shows a bend if this 
is the case for the kinetic temperature, but that the reverse is not true. 

6.2.2 Applications of the tests to models of nuclear fragmentation. 
Finite size scaling and microcanonical description 

The finite size scaling test for the determination of the order of a phase 
transition by means of the determination of critical exponents has been re- 
cently applied to a lattice model, the Ising Model with Fixed Magnetization 
(IMFM) ||197|| . Its canonical partition function reads 



Z v = ^2 ex P (P v o ^2 ai(r i) 5 m 



MM 
<n-n>i 



where Vq is a negative potential strength, M = 2p/po — 1 the "magnetization" 
related to the relative density p/po, M a fixed value of M , {<jj = ±1} related 
to the lattice site occupations {si = 0, 1} in the notation of subsection 5.3.5. 
The model is equivalent to an LGM except for the fact that the number 
of particles, A = £\ s iy is strictly fixed and, as a consequence, p/p is fixed 
if the volume V is fixed. It is a non trivial model whose properties have been 
investigated recently [j27^i ]. In the canonical ensemble the kinetic energy plays 



a trivial role, therefore there is no need to introduce it if one looks for critical 
properties of the system in this framework. 

The system governed by the IMFM shows a phase transition with a phase 
diagram (p, T) which has the same topology as LGM and related models, see 
f.i. Fig. [R]. However in the ordinary Ising model with zero magnetic field con- 
tribution, the magnetization jumps when the system enters the coexistence 
phase at T < T c . Since the present description keeps the magnetization fixed, 
the order of the transition is a priori not clear. 

The problem has been investigated numerically, first in the canonical 
framework for a system with fixed volume. The sensible quantity which has 
been worked out is the specific heat 

C v = 3V(< E 2 > - < E > 2 ) 
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where V is the volume and E the energy of the system. At the separation 
line, one looks for the asymptotic behaviour of Cy in the form 

C v oc (J3 - (3j h )~ a 

0Th _ j"-i^ pTh _ (^Th^-i^ rjiTh ^ e j n g fa e thermodynamic transition temper- 
ature of the infinite system for fixed density p and a a real positive exponent 
which may help to characterize the universality class to which the transition 
belongs. 

Using finite size scaling arguments and writing 

C V (L) = A + BL a/u (38) 

/3 c (L)-/3 c (oo) = AL- 1 /" 

where A, B, A are real constants and L is the linear size of the system, it is 
possible to obtain v and a/u, hence a, corresponding to the thermodynamic 
limit |197|| . Typical results for the evolution of Cy with the size of the system 
by identification of the maximum of this quantity for fixed p are shown in 
Fig. |25| for two values of p/p = (M + l)/2. In this figure a/u is read 
from the slope of the curves which should behave like rigorously straight 
lines if the scaling hypothesis (|38D is working. The straight line behaviour is 
observed for L large enough (L > 28), although there are deviations at close 
inspection. The values of v and a for p/p = 0.3 and 0.5 could be compatible 
with the universality class of the Ising model although this is not absolutely 
established since the behaviour for L = 48 could still be transitory. This 
value of L in 3d space requires already a large numerical effort. 

The non completely conclusive situation led to further investigations 



| 280|| . First, the topology of the system on the critical line is very simi- 
lar for different values of the density, showing a quite homogeneous pattern. 
Second, the evolution of the caloric curve shows a change of slope which gets 
more accentuated with increasing size of the system, with the presence of an 
inflexion point which is characteristic of a continuous transition. 

A sharper test consists of the use of the microcanonical ensemble. Here 
the kinetic energy part is taken into account and the procedure is described in 
subsection 5.3.6. As discussed above, a first order transition should manifest 
itself by a temperature which is a multivaluated function of energy, and 
this is not the case if the transition is continuous. This phenomenon was 
numerically checked on the Potts model which can effectively be of first or 
second order for different numbers of spin directions on a given site [ 1601 , 281 . 
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The application of the microcanonical algorithm to the IMFM does not show 
the expected "S"' curve as it would be the case for a first order transition, 
see Fig. ^6|. This fact does however not lead to a definite answer because it 
is expected that for a fixed volume the caloric curve in the infinite system 
never shows a plateau in the presence of a first order liquid-gas type phase 
transition but only a more or less steep change of slope, the plateau arises for 



fixed pressure as recalled in ref. [282]. Furthermore, as already mentioned, 



the finite system announces a transition which may not necessarily be the 
one which is really observed in the thermodynamic limit. 

Two conclusions can be drawn from this study. First, the present model 
deals with two constraints, fixed volume and fixed number of particles. These 
constraints may play a critical role in finite systems. This can be seen on the 
Potts model. For q > 5 where q is the number of spin directions on a site, 
the caloric curve should show an "S"' curve. When constraints are applied 
on the spin population of different spin directions at the different sites of the 
system the bend disappears ||280|| . 

The importance of constraints has of course also consequences for the 
interpretation of experimental results which depend of the physical conditions 
under which data are collected. Fixed volume is an illustration of the effect 
of constraints. 

6.2.3 Fluctuations and criticality 



As already seen above [|273|| , critical phenomena are strongly correlated 



with large fluctuations in the thermodynamic variables. This property has 
been exploited recently by Chomaz and Gulminelli |272| , |283| , [284|| in the frame- 



work of the microcanonical ensemble. These authors have explicitly shown 
on a simple model how the specific heat which can be defined in terms of 
energy fluctuations can reach negative values. 

A closed system with energy E t is divided into two subsystems, 1 and 
2, with energy E\ and E?, E t = E\ + E 2 - Then the energy distribution in 
subsystem 1 can be written as 



W(E t ) 
= Xexp[S 1 {E 1 )+S 2 {E t -E 1 )-S{E t )] (39) 

where N is a normalization factor, Wi(Ei) is the density of states of the 
system % at the energy Ei and Si(Ei) its entropy. W(E t ) is a convolution of 
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partial densities 



W(E t ) = / dE 1 exp[S 1 (E 1 ) + S 2 (E t - E 1 )) 
Jo 



At equilibrium W(E t ) is maximum for 



as. 



dE 1 



dS, 



E x 



dEo 



(40) 



Et-Ex 



where E\ corresponds to the stationary point where the equality is realized. 
Relation ( f40"D shows that the temperatures T% and T 2 are the same in both 
subsystems, T\ = T 2 = T as expected for a system in thermodynamic equi- 
librium. 

Close to the maximum of W(E t ) it is reasonable to parametrize P 1 * by 
a gaussian 



_2x-l/2 



Pf< (E ± ) = (27Taf)- l ^exp[-(E 1 -E i y/2a 



Then 



d 2 p^ t {E 1 



dEl 
From ©, for N = (27HX 2 )- 1 / 2 

9 2pBt 



"(27T) 



-1/2 



<J, 



(41) 



(42) 



dE\ 

where d 



21 Citf 2 



T 



<9 2 S, 
P 2 1 d E? 



(43) 



-i 



dP, 
dE t 



is the specific heat of the system i. Equating (|42| ) and ( f43|) and introducing 

C = Cx + C 2 (44) 

leads to 

C 2 



C 



Ci - o-f/T 2 
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If now subsystem 1 corresponds to the kinetic energy part of the sys- 
tem with energy E\. and specific heat C& and subsystem 2 to the potential 



contribution E p then 








C = 


ci 




Ck - crl/T 2 



(45) 

where C& = \A and A is the number of particles. The knowledge of the 
temperature T of the system and the kinetic energy fluctuations character- 
ized by a\ allow then to determine the heat capacity of the system in the 
framework of the microcanonical ensemble. It is clear that if <y\ gets larger 
than CkT , C gets negative which corresponds to an "S"' bend in the caloric 
curve. The analysis requires however some care. There is a priori no justifi- 



cation for the distribution (41) to be gaussian nor for ([0]) to be strictly true 
in the vicinity of a critical point. The lack of concavity of the entropy is the 
sign for it. There are however signs that, even if ( |I ) is not strictly correct, 



the variance of the real distribution can be assimilated to the variance of a 



gaussian distribution if the number of particles is not too small ||285|| . The 
correction to the heat capacity due to higher moment contributions to (|42D 
has been discussed in ref . | |272| . The kinetic energy distribution is in principle 
experimentally accessible. 



The relation (f44|) has been applied to lattice models of the same type as 



the one whose Hamiltonian is used in subsection 5.2.2 | 284 |. The partition 
function reads 



Z X (E) = Y,W V (E) 



e xv 



where Wy{E) is the density of states at energy E. The volume of the system 
is not fixed here but fluctuates. This is described by means of the introduction 
of a Lagrange multiplier A which is interpreted as a quantity proportional to 
a pressure P = \T\ where T\ = din Z\/dE. 

The equation of state for this model has been studied numerically in the 
3-dimensional space (T, E, A). One observes a backbending effect if one repre- 
sents P = XT as a function of the average volume < V > which is interpreted 
here as being the consequence of mass and energy conservation constraints. 
An "5"' bend appears also in the caloric curve on the (E, T) surface. For 
strictly fixed volume the caloric curve does not show any backbending. This 



is illustrated in Fig. |27| in which the specific heat is represented as a function 
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of energy using (^3). With this interpretation the puzzle concerning the order 
of the transition which was discussed in subsection 6.2.2 ||280|| is apparently 
explained, the shape of the caloric curve depends on the thermodynamic 
variables which are considered ||282|| . This message is of some importance in 
the framework of experimental analyses. 



6.3 Experimental construction of the caloric curve and 
determination of the specific heat 

Different attempts have been made since the pioneering work already 
mentioned above [ALADIN, INDRA, EOS, ... ] in order to construct the 
caloric curve starting from experimental data. 

The reaction 28 Si on 100 Mo leading to an incomplete fused system has 
been used in order to relate the excitation energy per particle e* of the system 
to a so called apparent temperature T app which is obtained from the slopes 
of the energy distributions of emitted light particles [ [286| . 



da 



dE, 



oc (E kin - B) exp 



kin 



(E, 



kin 



B) 



T, 



"/'/' 



where Ekin is the kinetic energy and B the Coulomb barrier. T app which is not 
the thermodynamic temperature T t h is shown to be a sensible quantity, both 
T app and T t h show the same trend as a function of the energy. The outcome of 
the experience is in reasonable agreement with calculations performed in the 
framework of the MMMC model. Since this model indicates the existence of a 
phase transition, it is tempting to conclude that it is also experimentally true 
since model calculations and experimental data have been worked out under 
the same physical conditions. It should nevertheless be mentioned that the 
results do not constitute a full proof because there remain some differences 
between the calculations and the data analysis. Former data from the Texas 



A & M-group [[260|1 are in agreement with those discussed here. 

The problem concerning the determination of the temperature has also 
been raised in the framework of a sequential decay model in which fragmen- 
tation evolves in time by means of successive binary break-ups of excited 
fragments [|287|| in a way which is similar to the approach of ref. [|94 |. The 
analysis [ p.42| , |288| of different definitions of the temperature leads to differ- 
ent results, due to the time-dependence of the process. It comes out that 
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the experimental caloric curve | 243 | can be reproduced over a large range of 
excitation energies. 

Very recently an attempt to detect the signal of a phase transition through 
the determination of the specific heat has been made by the Multics-Miniball 
collaboration at the cyclotron of the Michigan State University |289| . |290 



The excited system has been generated by means of Au on Au collisions at 
35 MeV-A. The angular distributions of fragments show isotropy which is a 



sign for possible equilibration |289|| . The conditions are such that the number 
of involved particles and the total energy are fixed, so that the arguments 
which characterize microcanonical ensembles can be applied. The correlation 
between the excitation energy E* and the temperature is fixed through the 
balance in energy at the freeze out 

M 

m + E* = ^( mi + ET l ) + E clb + ~(M - 1)T 

where mo is the initial mass excess, mi the mass excesses of the fragments, 
E* the excitation energy of the source, Ef 11 the internal excitation energy 
of fragment i, Ecib the Coulomb repulsion energy, M the multiplicity and T 
the temperature which is fixed by the balance equation. No collective energy 
is taken into account. A backtracing procedure is used in order to correlate 
the asymptotically observed events to the freeze-out configurations. The 
comparison of E*(T) with SMM predictions are in good though not perfect 
agreement. The heat capacity shows a sizable enhancement for temperatures 
lying between 4 and 6 MeV. 

In a further step the experiment has been analysed following the mi- 
crocanonical prescriptions of Chomaz and Gulminelli, see subsection 6.2.3 
||283| , |284|| . The total energy E t has been decomposed into a kinetic and a 
potential contribution 

Et = m + E* = E x + E pot 

M 

where E pot = ^m; + E C w 

i=l 

and Ecw is the Coulomb energy acting between the fragments at freeze-out. 
This energy has been evaluated by means of simulations in which primary 
fragments are randomly distributed in a freeze-out volume equal to three 
times the normal volume. Two extreme hypotheses have been used. Either 
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(1) all charged species including all light particles were considered as primary 
or (2) the totality of neutrons and light particles was distributed over the final 
fragments. These choices correspond to the procedures used in the framework 
of the SMM and MMMC models. In this way, one fixes the ensemble averaged 
quantities over sets of events < E\ >, < E\ >, <j\ and C\ — | < (M — 1) > 
where 

M 

< E, >= (J2^)t 2 + (-(M ~ 1))t 



=i 



and di is the level density of fragments i. 

The results corresponding to cases (1) and (2) are shown in Fig. |28. 
In both cases one observes the signs for a divergent behaviour of the heat 
capacity at two values of the energy and the expected negative contributions 
to the specific heat in the microcanonical framework. The stability of this 
result with respect to variations in the size of excitation energy bins, the 
assumptions on E*, E t and the choice of the a^'s have been checked. 

The same type of analysis has been performed with the INDRA detector 
on the outcome of the reaction Xe on Sn at 32 and 50MeV-A central col- 
lisions ||147| , [291f . In order to extract the kinetic energy fluctuations which 
enter the determination of the specific heat, it has been necessary to sub- 
stract the collective effects which are present and to take care of the particle 
emission which takes place during the expansion phase of the process. One 
of the two divergences which develop in the heat capacity is observed as well 
as the onset of negative values for this quantity. The second singularity is 
not seen because the excitation range does not cover the whole interval of 
energy corresponding to the coexistence of the phases. 
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7 Summary and conclusions 

In the present report we tried to summarize the present knowledge about 
the properties of excited nuclear matter which can be gained through the 
confrontation of simple classical models describing finite systems and exper- 
imental facts obtained from energetic collisions between finite nuclei. 

As an introduction we developed in section 2 some historical facts con- 
cerning the beginning of the activities in the field of nuclear fragmentation, 
both on the side of the experiment with the observation of power laws in the 
fragment size distributions, the introduction of the droplet model and the 
theoretical developments based on mean field approaches to nuclear matter 
and microscopic transport descriptions in the mean field approximation and 
beyond. 

We recalled that there are two main difficulties in the interpretation of 
experimental data obtained by means of energetic nuclear collisions. The 
first concerns the finite size of nuclei and the second the fact that the collid- 
ing objects evolve with time, from an initial compact system to an expanded 
system at the nuclear freeze-out where the nuclear interaction between frag- 
ments does no longer act, and finally to the asymptotic stage where the 
remnants of the reaction are identified by detectors. The second point leads 
to the question of the setting up of chemical and thermodynamic equilibrium. 
At lower energies there are signs which show that nuclear systems are able 
to absorb and spread sizable amounts of energy more or less uniformely over 
phase space. 

Thermodynamic equilibration is discussed in section 3 in which some of 
the efforts which have been made in order to get information about this point 
are developed. For the moment, equilibration seems to be a sensible concept 
since experimental results can be consistently interpreted in this framework, 
even though no strong and clear-cut proof for it exists up tp now. 

One of the first and more spectacular successes in the field concerns the 
multiplicities of fragments of different sizes which are amazingly well repro- 
duced by means of simple minded percolation models. This point is presented 
and discussed at length in section 4. Universality properties are clearly iden- 
tified and the intriguing power law behaviour which was observed and was 
of central interest in the early days appears again in a new context. Since 
percolation concepts are minimal information concepts it may appear that 
there is not much information to be gained from observables like fragment 
size distributions. The reasons for this may be related to the fact that highly 
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excited nuclei behave like disordered systems, and as it is the case for many 
other systems of this type in other fields of physics, the details of their dy- 
namical properties are of very restricted importance. Universality properties 
of physical systems are remarkable because they work as a link between ob- 
jects and phenomena which are seemingly very different from each other and 
at the same time are insensitive to the detailed information about the system. 

The success of percolation drew the attention towards other simple clas- 
sical models which have a strong kinship with percolation but are defined in 
terms of Hamiltonians. Some of them, so called lattice and cellular models, 
are described and analysed in section 5. Their simplicity and at the same 
time their universal character leads to the observation that they work as good 
tools for the investigation of the properties of fragmenting nuclei. Most of 
them are able to reproduce the properties related to fragment multiplicities. 
They show very useful for the investigations of effects related to finiteness in 
search for and analysis of phase transitions. They predict the possible exis- 
tence of two types of transitions which are not necessarily correlated in an 
evident way, a thermodynamic transition and a continuous percolation type 
transition which for systems with a density between nuclear density and half 
nuclear density occurs at temperatures which are higher than those which 
correspond to the thermodynamic transition. This point is not yet clearly un- 
derstood, although it does not lead, a priori, to any contradiction. It opens 
up the old question concerning the link with Fisher's droplet model. Nu- 
merical simulations on different models lead to very similar thermodynamic 
properties. 

The final part, section 6, deals more specifically with the problem of 
the existence and characterization of phase transitions and signs which show 
their existence in finite systems. For such types of analyses specific tools are 
needed. Some of them have been proposed very recently, at a time where sim- 
ilar considerations are developed in other fields of physics dealing with small 
systems. Some of these tools are explicitly described. Theoretical methods 
relying on a microcanonical description and aimed to detect transitions which 
have been developed in the framework of lattice models have been presented. 
Their application in the analysis of experimental data is described at the end 
of the section. 

A great wealth of information has been accumulated about the character- 
istics and behaviour of highly excited nuclei, their fragmentation properties 
and the existence of a genuine phase transition of a type which may be sim- 
ilar to the classical liquid-gas transition. The last point remains however to 

103 



be confirmed. Many pieces of the puzzle have been brought together, how- 
ever the whole picture is not yet totally clear. Experimental confirmation of 
the transition has to be brought, scaling properties further investigated, the 
question of the possible coexistence of two types of transitions elucidated. 
Quantum mechanical aspects which have been quoted and the importance 
of quantum effects have to be further worked out and analysed. Last, there 
remains the difficult problem of the possible off-equilibrium character of the 
fragmentation process of a finite critical system. The study of critical off- 
equilibrium dynamics is a relatively recent field of research, still in its in- 
fancy. It is a field to which the nuclear physics community interested in the 
behaviour of highly excited and decaying systems of particles could certainly 
bring its contribution. 
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FIGURES 
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Figure 1: Equation of state relating the pressure to the density in nuclear 



matter |25|1 . The curves represent isotherms. The spinodal region is indicated 
by dashed curves. 
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Figure 2: Mass distributions for an expanding system of A = 100 particles 
for different initial temperatures | |S4j| . The dots are obtained by means of 
CMD calculations, the lines are fits using Fisher's formula. 
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Figure 3: Histogram of mass yields obtained by means of the rate equations 
(0) after t = 4 x 10 _23 s [79]. The dots are experimental results of the reaction 
p + Kr at E = 80 - 350 Gel/. 
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successive break-ups in the fragmentation process with increasing excitation 



energy 104 



125 



A 
V 




I I I I I I I I I I I I I I I I I 

1000 hO X-LCP, • IMF,60AMeV 

" D X-LCP, ■ IMF,50AMeV 

e _ rt A X-LCP, A IMF,40AMeV 

O 800 :* X-LCP, • IMF,JOAMeV 

600 h 







.oc 



x 



_L 



_L 



200 400 



300 1Q00 1200 1400 1600 



Et (MeV) 



Figure 5: Top panel : average number of IMFs (solid symbols) and LCPs 
(open symbols) as a function of the total transverse energy E t . Bottom 
panel : corresponding average transverse energies as a function of E t . The 
different symbols correspond to different bombarding energies indicated in 
the bottom panel |p.43|j . See text. 
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Figure 6: Angular distribution of light particles (Z = 1,2) in the flow angle 
6 > 60° corresponding to events generated by a unique source in the reaction 



Xe + Snat 50MeV-A |47 |. 
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Figure 7: Plot of the (a, s) plane where the degeneracy g(a, s) is non zero. 
For fixed A (= A\ or A 2 ) only the hatched areas contribute to the mass 
spectrum ||168|| . 
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Figure 8: Mass yields obtained with the parametrization of the bond proba- 
bility p(b) given in the text for (a) p + Ag at 11.5 GeV, po = 0.65; (b) p + Ta 
at 5.7 GeV, p = 0.69; (c) p + 197 Au at 11.5 GeV, p = 0.70 pTlj . 
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Figure 9: In S3 as a function of MS2 for single events. Part (a) corresponds 
to the break-up of gold nuclei, part (b) to a cubic bond simulation with 216 
sites and bond probabilities < p < 1, [|174|| . 
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Figure 10: The quantity 72 defined in (£2j). The upper part corresponds to 
the fragmentation of gold and the lower part to percolation simulations, ||1 75| . 
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Figure 11: Size of the largest fragment for different events as a function of 
£2. (a) Experimental results obtained from 376 events, (b) Bond percolation 
results in a cubic volume with 216 sites, 4000 realizations ||174j . 
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Figure 12: Left: measured mean multiplicities < M IMF > as a function 
of Z bound for 238 U on 197 Au (circles), 197 Au on 197 Au (squares) and 129 Xe on 
197 Au (triangles) for E = 600 MeV-A. Right: same quantities but rescaled by 
normalization with respect to the number of charges Z p of the projectile ||1 79| . 
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Figure 13: Average values of different observables defined in the text as a 
function of Z bound ||180| . 
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Figure 14: Phase diagram showing a first order transition for eight fluids 
indicated on the figure |192| . The collapse of the coexistence curve shows the 
universality character of the phenomenon. 
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Figure 15: Multiplicities of clusters made of parallel spins in the Ising model 
for (1) short range and (2) long range interactions and different temperatures 
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Figure 16: Typical phase diagram corresponding to a lattice gas model |20£ 
See discussion in the text. 
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Figure 17: Exact (full lines) and mean field (dotted lines) chemical potentials 
as a function of pressure for different temperatures |223|| . 
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Figure 18: Comparison of percolation results (upper part) with A2 algorithm 
(lower part). |223j. 
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Figure 19: Different observables represented as a function of Zb oun d. The ex- 
perimental average largest cluster < Z max >, IMF multiplicities < Mimf > 
and He multiplicities < M# e > represented by dots are compared to calcu- 
lations performed with the AO algorithm, see text [224]. 
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Figure 20: Mass distributions for a 6 3 system. The upper curve corresponds 
to a calculation using the AO algorithm, the lower corresponds to the simu- 



lation of collective effects, see text [224]. 
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Figure 21: Phase diagrams for different types of interactions, see explanations 



in the text 235 
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Figure 22: Caloric curve and specific heat for a finite size system, L = 8 3 , and 



density p/po = 0.5. Canonical and microcanonical simulations coincide ||2 1 1 



143 



10 



a - 



LU 



6 - 






+ - 



CELLULAR MODEL L=S 
No Coulomb interaction 



T£thermo) 
T(droplet) 







y 



n ! i i i I i i i I i i i I i i i I i i l 

02 DA o,e 0,3 1 

DENSfTY 



Figure 23: The phase diagram obtained with the cellular model described in 
the text ||241|| . See discusion in the text. 
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Figure 24: Isotope temperature T HeLi (full triangles) and temperatures ob- 
tained by means of excited state population rate measurements (open sym- 
bols) as a function of the excitation energy per nucleon for 197 Au on 197 Au 
at E= 1000 MeV- A [g2|. 
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Figure 25: Fit of Cy ax as a function of the linear size L of the system 
for different densities. Crosses indicate different sizes. The largest size is 




Figure 26: Caloric curve calculated for a 3D IMFM for a system of linear 
size L = 10, p/po = 0.3 (upper part), p/po = 0.5 (lower part). Canonical 
and microcanonical calculations are indistinguishable. See comments in the 
text |280H . 
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Figure 27: Caloric curve, specific heat and kinetic energy fluctuations at 
constant pressure (left part) and constant volume (right part) calculated in 
the framework of the microcanonical ensemble |[284| . See comments in the 
text. 
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Figure 28: Experimental heat capacities for the two extreme hypotheses 
which where made in order to evaluate E C w (see text). The dotted curve 
corresponds to the Fermi-gas calculation 
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